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Abstract 

A  high  speed  durable  ion  probe  based  pressure  sensor  is  being  investigated  for  use  in 
pulse  detonation  engines.  The  environment  eneountered  in  such  engines  necessitates  high 
temperature  and  durable  (vibration  resistant)  devices.  Traditional  pressure  sensors  can  be 
used  however,  various  methods  and  materials  used  to  protect  the  sensors  dampen  and 
reduce  the  pressure  wave  allowing  for  qualitative  results  only.  An  alternative  transient 
pressure  sensing  method  is  investigated  for  pressures  behind  a  hydrocarbon  flame  in  the 
pulse  detonation  engine.  Hydrocarbon  flames  generate  ions  that  are  quenched  by 
collisions  with  other  species  and  walls.  As  the  collision  rate  is  a  function  of  pressure,  so 
too  is  the  ion  decay  rate.  The  ion  decay  rate  is  measured  using  an  ion  probe  that  is  well 
suited  for  high  temperature  flow,  has  no  moving  parts,  and  is  inexpensive.  Similar 
systems  have  been  used  to  determine  multiple  combustion  conditions  in  automobile 
engines.  This  investigation  builds  upon  these  capabilities  to  examine  the  quantitative 
pressures.  The  ion  probe  measures  the  ionization  in  the  form  of  a  small  current.  The 
strength  of  the  ion  current  indicates  the  strength  of  the  ionized  field  which  decays 
according  to  pressure.  An  experiment  was  devised  to  correlate  the  ion  current  decay  rate 
with  the  pressure.  A  correlation  has  been  established  showing  pressure  is  a  function  of 
the  ion  current  decay  rate.  This  investigation  shows  a  viable  alternative  method  for 
measuring  pressure  in  the  pulse  detonation  engines  although  additional  work  is  required 
to  improve  the  accuracy  of  the  method. 
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ION  BASED  PRESSURE  SENSOR  FOR  PULSE  DETONATION  ENGINES 


I,  Introduction  and  Overview 


1.1  Chapter  Overview 

This  chapter  describes  the  need  for  an  improved  pressure  sensor  for  the  pulse 
detonation  engine.  The  methodology  for  making  this  improvement  is  also  explained. 

1.2  Motivation 

The  United  States  Air  Force,  along  with  other  organizations,  is  currently 
investigating  the  Pulse  Detonation  Engine  (PDE)  as  a  future  propulsion  system.  While 
not  a  new  concept,  the  engine  is  still  in  developmental  stages.  The  PDE  theoretically 
offers  higher  efficiency  with  less  complexity  and  lower  weight  than  the  turbofan  engines 
in  use  today.  In  addition  to  air-breathing  cycle,  the  PDE  can  also  operate  as  a  rocket  cycle 
termed  the  Pulse  Detonation  Rocket  Engine  (PDRE).  The  PDE  and  PDRE  are  also 
attractive  due  to  their  large  flight  envelope:  from  static  up  to  around  Mach  Number  5. 
Turbofan  designs  are  typically  limited  to  Mach  Number  2  or  3.  On  the  other  hand, 
ramjets  and  scramjets  require  supersonic  speeds  in  order  to  start  producing  thrust.  The 
large  flight  envelope  of  the  PDE  eliminates  the  need  for  any  boosters. 

The  basis  for  the  PDE  is  the  higher  efficiency  of  a  detonation  combustion  process 
compared  to  the  constant  pressure  deflagration  process  used  in  conventional 
turbomachinery  based  air-breathing  engines  of  today.  This  efficiency  comes  from  the 
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near  constant  volume  process  and  the  fact  that  the  PDE  does  not  require  the  working  fluid 
to  be  compressed  prior  to  heat  addition.  Although  unsteady,  this  process  closely  follows 
the  thermodynamics  of  a  Humphrey  constant-volume  cycle  [6].  Simple  theoretical 
calculations  show  the  efficiency  of  the  Brayton,  Humphrey,  and  Chapman- Jouguet 
Detonation  cycles  to  be  27%,  47%,  and  49%  respectively  [6].  Compared  to  the  constant 
pressure  Brayton  cycle,  the  Humphrey  cycle  achieves  higher  efficiency  by  creating 
higher  temperatures  at  lower  entropy.  In  addition  to  the  clear  thermodynamic  advantages, 
the  PDE  also  has  the  potential  to  reduce  cost  and  enhance  performance  without  the  heavy 
turbomachinery  in  conventional  air-breathing  engines. 

Conventional  turbomachinery  based  engines  use  a  steady  process  of  compression, 
heat  addition,  and  expansion  to  generate  thrust.  The  pulse  detonation  engine  generates 
thrust  through  an  entirely  different  unsteady  process.  PDE  is  similar  in  many  ways  to 
internal  combustion  (IC)  engines.  Eike  the  IC  engine,  the  PDE  fills  a  tube  with  air  and 
then  adds  fuel  creating  a  near  stoichiometric  mixture.  In  the  IC  engine,  the  piston 
compresses  the  mixture  and  initiates  deflagration  using  a  spark  plug.  In  the  PDE,  no 
compression  is  required.  The  fluid  can  also  be  ignited  with  a  spark  plug,  but  deflagration 
instead  transitions  to  detonation  as  the  combustion  moves  down  the  tube.  This  detonation 
wave  is  the  basis  of  the  PDE.  Detonation  by  nature  is  an  unsteady  process  where  the 
wave,  according  to  Chapman-Jouguet  theory,  travels  at  supersonic  speeds  relative  to  the 
unbumed  fuel-air  mixture.  The  PDE  takes  advantage  of  this  unsteady  process  by 
employing  multiple  detonation  tubes  similar  to  multiple  cylinders  in  an  IC  engine.  Each 
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tube,  similar  to  the  4  stroke  eycle  in  IC  engines,  is  either  being  filled  with  the  fuel-air 
mixture,  detonating  the  mixture,  blowing  down,  or  purging  the  exhaust. 

In  order  for  a  sueeessful  detonation  near  Chapman- Jouguet  predicted  speeds,  the 
combustion  must  produce  a  strong  shock  wave  that  travels  down  the  tube.  This  shock 
wave  increases  the  temperature  and  pressure  of  the  fuel-air  mixture.  After  a  short 
induction  period,  this  mixture  combusts  or  detonates  in  a  thin  region  behind  the  shock 
wave.  The  detonation  then  creates  the  even  higher  temperatures  and  pressures  needed  to 
sustain  the  shock  wave. 

Each  detonation  wave  produces  a  small  amount  of  thrust  based  on  the  diameter  of 
the  tubes,  the  speed  of  the  detonation  wave,  and  the  pressure  behind  the  wave.  Since  the 
PDE  is  unsteady,  the  rate  of  firing  each  tube  directly  affects  the  generation  of  thrust. 
Substantial  gains  in  thrust  can  be  achieved  by  increasing  the  cycle  frequency.  At  higher 
frequencies,  however,  timing  becomes  critical  to  successful  detonations.  Although  the 
deflagration  to  detonation  transition  has  been  heavily  researched,  in  practice  wave  speeds 
near  Chapman- Jouguet  theory  are  not  always  realized.  Often,  weaker  shock  waves  are 
formed  resulting  in  substantially  slower  wave  speeds.  These  weak  detonations  greatly 
reduce  thrust  of  the  engine. 

Despite  recent  progress,  significant  challenges  remain  before  reliable  PDE 
operation  with  practical  fuels  is  realized  [11].  Eurther,  the  cycle  creates  higher 
temperatures  than  the  Brayton  cycle  leading  to  high  heat  loading  [10].  As  PDEs  increase 
their  cycle  frequency,  heat  loads  increase  [10]  thus  heat  related  problems  will  only 
worsen.  These  high  temperatures  limit  the  diagnostic  tools  available  to  researchers. 
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Specifically,  conventional  piezoelectric  based  pressure  transducers  are  ill-suited  for  the 
high  temperatures  and  harsh  vibratory  environment  within  the  PDE.  A  variety  of 
techniques  can  be  used  to  increase  the  useful  limits  of  the  piezoelectric  pressure 
transducers.  Each  of  these  techniques  have  disadvantages  that  often  skew  the  results.  Eor 
example,  protective  ablative  coatings  on  the  pressure  transducers  improve  the  resistance 
to  the  harsh  environment  but  reduce  the  sensitivity.  These  coatings  reduce  the 
effectiveness  of  the  pressure  transducer  as  a  quantitative  instrument  because  of  the 
inherent  dampening  of  the  materials.  Accurate  compensation  for  the  dampening  effects  is 
not  feasible  due  to  the  variability  in  the  thickness  of  the  material  as  well  as  the  ablation 
rate  itself.  Eor  single  firings  of  a  detonation  tube  this  ablation  can  be  measured  and  added 
as  a  correction  factor  to  the  pressure  measurement.  In  steady  operation,  measuring  the 
ablation  of  the  protective  material  is  not  feasible.  A  durable,  quick  response  quantitative 
pressure  sensor  is  needed  to  optimize  the  PDE  during  development,  and  also  to  provide 
feedback  for  engine  control. 

1,3  Method 

IC  engines  have  also  encountered  similar  problems  with  measuring  pressure 
within  the  cylinder.  Sensing  the  cylinder  pressure  enables  tighter  control  of  equivalence 
ratio  (cp)  leading  to  reduced  hydrocarbon  emissions  [15].  Equivalence  ratio  by  definition 
is  the  air  to  fuel  ratio  divided  by  the  air  to  fuel  ratio  at  stoichiometric  conditions.  While 
various  techniques  exist  for  developmental  engines,  modifying  production  engines  to 
include  a  reliable  pressure  sensor  is  not  practical.  Production  engines  have  no  place  to 
install  extra  high-cost  sensors.  Using  these  conventional  sensors  can  change  the 
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combustion  in  the  cylinder  and  require,  similar  to  the  PDE,  a  eomplex  eooling  system 
[21],  The  sensors  and  modifieations  are  not  praetieal  for  long-term  use  needed  for 
produetion  engines.  In  order  to  solve  this  problem,  several  teehniques  have  been 
developed  to  extraet  information  from  the  ionization  resulting  from  the  combustion 
process.  By  plaeing  a  small  direct  current  (DC)  voltage  aeross  the  gap  after  diseharge,  the 
spark  plug  acts  as  an  ion  sensor.  This  ionization  then  produees  a  small  ion  eurrent  aeross 
the  spark  plug  gap  after  the  ignition.  Without  any  modifications  within  the  eylinder, 
information  ean  be  extraeted  from  the  ion  eurrent.  To  date,  the  spark  plug  has  eontrolled 
the  equivalence  ratio  [15],  deteeted  misfire  [3],  and  eontrolled  knock  [3].  Additionally, 
the  spark  plug  has  been  employed  to  measure  pressure  in  the  cylinder  [19].  In  short, 
internal  combustion  engines  have  utilized  the  ion  current  across  the  spark  plug  to 
measure  several  important  eonditions. 

Applied  to  the  PDE,  the  spark  plug  already  aets  as  a  rugged  ion  sensor  to  measure 
wave  speed  [22,  23].  Two  spark  plugs  with  an  applied  DC  voltage  are  inserted  a  known 
distanee  apart  in  the  PDE  tube.  The  measurement  is  made  by  determining  the  time  delay 
between  the  voltage  diseharges  of  two  spark,  resulting  in  a  simple,  but  highly  useful 
method  for  determining  average  detonation  wave  speed  [22,  23]. 

The  spark  plug  has  already  proven  its  durability  to  the  harsh  PDE  environment. 
Extending  the  use  of  the  spark  plug  to  measure  pressure  in  the  PDE  is  a  logieal 
improvement.  By  utilizing  the  advancements  in  IC  engines,  the  spark  plug  ean  be 
employed  to  a  mueh  greater  extent  in  the  PDE  and  become  an  additional  pressure 


5 


transducer.  Before  this  extension  can  be  put  to  use,  the  underlying  theory  must  be  tested 
for  conditions  in  the  PDE  and  several  engineering  ehallenges  overcome. 

1,4  Thesis  Content 

This  thesis  covers  the  experimental  work  completed  in  adapting  the  spark  plug  as 
a  pressure  sensor  in  the  PDE.  Previous  work  and  the  theoretical  basis  of  this  investigation 
are  deseribed  in  Chapter  II.  The  progress  made  in  using  the  spark  plug  as  an  ion  sensor  in 
internal  eombustion  engines  is  leveraged  and  applied  to  this  investigation.  Based  on  the 
previous  work,  predictions  are  made  on  how  the  deeay  rate  of  the  ion  current  is  a  function 
of  pressure.  In  order  to  test  this  prediction,  an  experimental  approaeh  is  devised.  The 
details  of  the  approach  including  the  instrumentation  and  data  eapture  are  described  in 
Chapter  III.  The  raw  results  from  this  testing  are  shown  in  Chapter  IV  along  with  a  short 
discussion  on  the  observed  phenomena.  The  data  is  then  analyzed  and  redueed  in  Chapter 
V.  An  error  analysis  is  accomplished  to  determine  the  accuraey  of  the  data.  The 
predictions  of  Chapter  II  are  compared  to  the  analyzed  data  in  Chapter  V  and  the  overall 
accuracy  of  this  method  is  addressed.  Conelusions  of  this  investigation  are  described  in 
Chapter  VI. 
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II,  Background  and  Theory 


11.1  Chapter  Overview 

This  chapter  describes  the  work  of  others  in  understanding  the  theory  on  how  to 
determine  pressure  from  the  ion  sensor.  This  knowledge  is  then  applied  to  the  problem  at 
hand  by  predicting  how  the  sensor  will  fimotion. 

11.2  Ion  Formation 

It  is  well  known  that  hydroearbon  flames  have  eonductive  properties. 
Considerable  researeh  over  several  deeades  has  investigated  the  formation  of  ions  in 
flames  [8].  The  formation  of  these  ions  within  the  flame  is  attributed  to  the  ehemi- 
ionization  reaetions  [4].  These  reaetion  are  not  initial  reaetants  and  produets  but 
intermediate  short  lived  species  of  the  eombustion  proeess.  One  important  example 
being  [5]; 


CH  +  O  ^  CHO^  +  e' 


(1) 


Other  typieal  reaetions  include  [3,21]; 


CHO^  +H2O  ^ 

"  CO  +  H30^  +  e 

(2) 

CH  +C2H2  - 

CsHb"^  +  e' 

(3) 
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Many  other  reactions  are  important  in  chemi-ionization  and  models  for  ionization 
can  be  very  complex.  Because  these  reactions  and  species  can  vary  greatly  depending  on 
the  initial  and  local  conditions,  the  chemi-ionization  process  can  be  extremely  complex. 
Extensive  research  has  investigated  these  species  and  reactions.  In  a  laboratory  setting 
important  reactions  and  species  can  be  identified.  Unknowns  and  local  variability  in  the 
PDE,  however,  preclude  the  detailed  examination  of  individual  species  and  reactions. 
Instead,  the  ions  will  be  considered  at  a  global  level  for  the  development  of  a  useful 
pressure  sensor  in  the  PDE. 

Under  local  thermal  equilibrium  conditions,  the  ion  concentration,  as  a  function  of 
temperature,  is  given  by  Saha’s  equation: 


=  2 


~1? 


B.. 


B 


exp 


i-\ 


kT 


(4) 


The  reactions,  as  shown  by  the  example  in  Eq.  (1),  in  a  hydrocarbon  flame,  place  the  ion 
concentration  at  super-equilibrium  levels. 

11,3  Ion  Decay 

The  net  rate  change  of  ion  concentration  is  the  difference  between  the  production 
and  recombination  rates  of  reaction.  Typical  recombination  reactions  are  given  by  [3,  5]: 


+  e'  ^  H2O  +  H 
^  OH  +  2H 


(5) 
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The  super-equilibrium  ion  coneentration  will  naturally  deeay  to  the  levels  governed  by 
Eq.  4  through  moleeular  eollisions.  As  the  reeombination  reaetion  has  a  moleeularity  of 
two,  the  rate  of  ion  deeay  is  therefore  dependent  on  the  square  of  pressure  [4], 

By  measuring  the  ion  density  as  a  funetion  of  time  it  is  possible  to  determine  the 
pressure.  While  the  pressure  eould  be  theoretieally  determined  in  any  region  of  ion 
eoneentrations  above  equilibrium,  praetieality  limits  the  regions  for  useful  measurement. 
While  ion  production  can  be  predicted  in  a  tightly  controlled  laboratory  setting  where 
highly  sensitive  initial  and  local  conditions  can  be  determined,  implementing  these 
measurements  into  the  PDE  is  not  practical.  The  pressure  measurement  can  be  simplified 
if  the  decay  is  observed  well  past  the  ion  generating  reaction  front  so  that  only  the 
recombination  reaction  rate  need  be  considered. 

11,4  Internal  Combustion  Engines 

The  various  methods  in  IC  engines  have  already  used  these  chemi-ionization 
relationships  to  successfully  measure  various  properties  of  combustion  [3,  17,  21].  Some 
methods  refer  to  a  flame  resistance  instead  of  an  ion  current.  The  flame  resistance  is 
inversely  proportional  to  the  ion  current  based  on  Ohm’s  Eaw: 

V 

R  =  j  (6) 


where  V  is  the  voltage,  I  is  the  current,  and  R  is  the  flame  resistance. 
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Since  the  IC  engine  methods  use  the  spark  plug  both  for  eombustion  ignition  and 
ionization  detection,  careful  eonsideration  must  be  given  to  prevent  the  ignition  system 
from  interfering  with  the  ion  sensing.  The  PDE  setup  ean  avoid  some  of  the 
eomplications  of  the  IC  systems  by  separating  the  ignition  and  ion  sensing  functions. 
Many  significant  results  ean  be  direetly  ineorporated  for  use  in  the  PDE. 

11.4.1  Ignition  Spark 

In  some  eases,  the  ion  eurrent  may  be  obseured  by  the  spark.  The  spark  is  the 
result  of  the  breakdown  of  the  loeal  air  and  fuel  mixture  into  plasma  by  a  strong  eleetrie 
field.  During  the  spark,  the  eleetrie  field  and  plasma  will  dominate  the  current 
measurement  resulting  in  large  current  variations  not  based  on  ion  produetion.  The  sensor 
values  during  the  spark  must  therefore  be  disearded  when  solely  trying  to  measure  the 
current  due  to  ionization.  A  short  time  period  after  the  spark,  the  measured  eurrent  ean  be 
eonsidered  dominated  by  the  ionization  of  the  mixture. 

11.4.2  Sensor  Configuration 

Some  results  from  IC  engines  ean  be  direetly  applied  to  the  PDE.  Applying  a 
positive  DC  voltage  eaptured  a  larger  quantity  of  ions  than  a  negative  DC  bias  beeause  of 
the  higher  mobility  of  the  electrons  compared  to  the  positive  ions  [21].  Eurther,  the 
deteetion  sensitivity  improved  when  inereasing  the  surfaee  area  of  the  eenter  eleetrode  on 
the  spark  plug  [17]. 
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11,4.3  Combustion  Conditions 


When  testing  several  Air/Fuel  ratios,  results,  for  an  engine  running  at  2570 
revolutions  per  minute  (rpm),  showed  the  flame  resistanee  was  at  a  minimum  when  the 
peak  internal  pressure  was  at  a  maximum  [17].  A  low  flame  resistanee  eorresponds  to  a 
high  eurrent  through  Ohm’s  Law  and  shows  that  aetual  ionization  levels  are  highest  near 
stoiehimetrie  eonditions.  At  this  engine  speed  the  eombustion  and  post-eombustion  zones 
oeeur  on  the  order  of  12  milliseeonds  (ms)  assuming  a  ehange  in  erank  angle  of  90 
degrees.  Results  also  showed  that  the  intake  pressure  does  not  ehange  the  flame 
resistanee  signifieantly  [17].  This  finding  reinforces  the  fact  that  the  ionization  is  due  to 
the  chemi-ionization  and  not  initial  pressure. 

Additional  work  examined  the  use  of  ionization  current  to  adjust  timing  in  an  IC 
engine  [7].  One  ionization  measurement  system  is  already  in  use  in  a  SAAB  engine  [7]. 

IC  engines  use  a  peak  pressure  algorithm  for  ignition  timing  [7].  These  algorithms  are 
constrained  by  the  thermal  and  high  pressure  limits  of  pressure  sensors  [7]  and  could  be 
improved  by  using  the  ion  current  across  the  spark  plug  as  a  feedback  sensor  to  determine 
peak  pressures.  Ionization  current  can  be  affected  by  temperature,  air-fuel  ratio,  time 
since  combustion,  exhaust  gas  recycling,  fuel  composition,  engine  load,  etc  [7].  Despite 
complications,  results  show  typical  ionization  curves  for  ignition,  flame  front,  and  post 
flame.  In  the  post  flame  region,  relatively  stable  ions  follow  the  cylinder  pressure  trend 
[7].  NO  was  found  to  be  a  contributor  to  the  post-combustion  ionization  because  of  the 
low  ionization  energy  [7].  A  Gaussian  function  for  the  ion  current  was  developed  based 
on  the  pressure  [7].  Problems  arise  when  trying  to  extract  pressure  information  from 
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ionization  current  [7].  A  peak  pressure  seareh  is  not  feasible  since  the  flame-front  often 
has  more  than  one  peak  and  the  post- flame  zone  doesn’t  have  a  peak  [7].  To  address  the 
problem  of  flame  fronts  with  two  peaks,  two  Gaussian  models  were  used  for  the  flame 
front  and  one  for  the  post  flame  phase  [7].  This  teehnique  captured  the  structure  of  the 
ionization  eurrent  although  quantitative  eomparison  was  not  provided  [7]. 

Relatively  few  experiments  have  been  conducted  in  a  combustion  bomb  whereas 
most  experiments  have  been  conducted  using  gasoline  IC  engines.  A  constant  volume 
eombustion  bomb  experiment  was  conducted  to  investigate  the  use  of  the  spark  plug  as  a 
eombustion  probe  mainly  to  estimate  eombustion  quality  [2],  The  approaeh  only 
addressed  the  question  of  combustion  quality  and  not  the  underlying  combustion 
phenomena  [2],  Several  signal  proeessing  methods  were  used  in  this  analysis.  Pattern 
recognition  and  classifier  design  were  used  to  perform  signal  classification  [2]. 
Classification  also  was  accomplished  by  artifieial  neural  networks  using  Matlab  [2].  In 
addition  to  the  combustion  bomb  experiment,  tests  were  also  condueted  with  internal 
eombustion  engines.  This  experiment  eoncludes  that  a  nonlinear  relationship  exists 
between  ionization  eurrent  and  combustion  quality  [2].  Results  showed  that  the  spark 
plug  can  be  a  reasonable  ionization  probe  as  long  as  efficient  signal  proeessing 
algorithms  are  used  [2]. 

Further  work  shows  that  the  ionization  elosely  follows  the  pressure  variation  in 
time  seales  of  milliseeonds  in  a  combustion  bomb  setup  [1].  Temperature  is  assumed  to 
be  a  known  function  of  time.  In  a  combustion  bomb,  the  rate  of  rise  of  temperature 
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determines  the  maximum  value  of  the  ionization  current  [1],  The  current  is  also  very 
sensitive  to  the  air-fuel  ratio.  At  high  air-fuel  ratios,  the  current  quickly  decreases. 

11,4,4  Work  of  Saitzkoff  et  al, 

Saitzkoff  et.  al.  [18]  investigated  the  use  of  the  spark  plug  as  an  ionization  sensor 
for  internal  combustion  engines.  Their  work  assumes  thermodynamic  equilibrium 
conditions  after  complete  combustion  where  the  gas  is  undergoing  adiabatic  expansion. 
The  test  V6  engine  ran  at  1300  rpm  at  full  load.  At  this  engine  speed  the  combustion  and 
post-combustion  zones  occur  on  the  order  of  23  ms  assuming  a  change  in  crank  angle  of 
90  degrees.  By  applying  a  80  volt  DC  to  the  spark  plug  the  measured  voltage  was 
converted  into  current  across  a  known  22  kiloOhm  (kD)  resistance  with  an  estimated 
error  of  5%  [18].  An  adiabatic  maximum  flame  temperature  of  2800  Kelvin  (K)  was 
assumed  at  a  maximum  pressure  of  5.7  MPa  (56.25  atm)  [18].  Saitzkoff  et  al.  assume 
Nitric  Oxide  (NO)  to  be  the  dominate  ionization  molecule  due  to  a  low  ionization  energy 
of  9.27  electron  Volts  (eV)  [18].  NO  is  formed  by  means  of  the  extended  Zeldovich 
mechanism  [14]. 


0  +  N2  ^  NO  +  N 

(V) 

N  +  O2  ^  NO  +  0 

(8) 

N  +  OH  ^  NO  +  H 

(9) 

NO  can  also  be  formed  by  the  low  temperature  “prompf  ’  or  Fenimore  NOx  mechanism. 
High  NO  formation  rates  exist  near  the  combustion  zone  due  to  super-equilibrium  levels 
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of  O  and  OH  radicals  [18].  Since  the  eombustion  zone  is  thin  at  these  high  pressures, 
Saitzkoff  et  al.  assume  the  formation  of  NO  near  combustion  to  be  small  compared  to 
formation  in  the  post-eombustion  zone  and  therefore  unimportant  in  their  model  [18]. 

The  NO  concentration  in  the  post-eombustion  zone  is  assumed  to  be  1%  [18].  Based  on 
this  assumption,  Saitzkoff  et  al.  claim  that  the  source  of  the  free  eleetrons  is  not  ehemieal 
reactions  but  thermal  ionization  [18]. 

Using  the  thermal  ionization  assumption,  Saitzkoff  et  al.  derived  an  ionization 
model  using  Saha’s  equation  (4),  the  ionization  ratio  of  the  particles,  and  the  electron 
drift  velocity  [18].  Within  this  model,  the  normalized  current  and  pressure  values  were 
related  by  the  following  [18]; 


1  3r 

2  4  y 


-exd 


P 

V  max  J 


ff 


E, 


2kT. 


f  \ 


r-1 


-1 


y^VV-^maxy  y 


(10) 


Calculated  pressures  were  correlated  to  the  post-flame  ionization  peak  due  to  NO 
production  [18].  A  low  signal-to-noise  ratio  required  the  current  to  be  filtered  before 
making  caleulations  and  data  was  averaged  over  50  cyeles  [18].  Although  equation  (6)  is 
sensitive  to  the  temperature  because  of  the  exponential  term,  Saitzkoff  et  al.  found  the 
experimental  relative  values  to  be  slightly  higher  than  predicted  but  still  in  fair  agreement 
[18]. 
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Saitzkoff  et  al.  state  their  ionization  model  is  only  applicable  to  the  post¬ 
combustion  zone  [18]  but  contradict  themselves  by  applying  the  model  to  the  full  crank 
angle  limits.  This  limitation  causes  poor  prediction  of  the  ion  current  during  combustion 
as  expected  but  still  agrees  with  the  post-flame  ionization  peak.  At  these  high  load  test 
cases,  high  temperature  creates  the  needed  activation  energy  for  the  NO.  The  low  engine 
rpm  allows  NO  the  relatively  long  chemically  kinetic  formation  time.  The  high  pressure 
produces  a  high  partial  pressure  of  NO  thus  increasing  the  ion  density.  At  these  specific 
conditions  the  assumption  that  post-flame  ionization  is  dominate  is  valid.  Despite 
neglecting  chemi-ionization,  an  additional  smaller,  although  still  prominent,  ionization 
peak  occurs  before  the  post-combustion  peak  [18].  At  lower  engine  loads  the  lower 
temperatures  will  decrease  the  level  of  thermal  ionization  and  could  allow  the  chemi- 
ionization  peak  within  the  combustion  zone  to  become  dominate.  Depending  on 
combustion  conditions  either  or  both  ionization  peaks  may  be  important. 

11,4,5  Improved  Model  of  Saitzkoff  et  al. 

In  follow-on  work,  Saitzkoff  et  al.  [19]  sought  to  improve  the  results  of  trying  to 
predict  pressure  with  the  ion  current  across  the  spark  plug.  Saitzkoff  et  al.  [19]  again 
focused  on  the  post-flame  zone  where  the  gas  species  are  assumed  to  be  in  chemical 
equilibrium  and  only  thermodynamic  conditions  are  changing.  They  identified,  however, 
an  additional  smaller  ionization  peak  above  thermal  ionization  levels  due  to  the  chemi- 
ionization  processes  at  the  flame  front  [19].  They  relax  previous  assumptions  to  allow 
species  with  low  ionization  energies,  such  as  long  lived  hydrocarbons,  to  chemically  react 
and  influence  the  ion  current  in  addition  to  thermal  equilibrium  levels  [19].  A  zero- 
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dimensional  chemical  kinetic  model  was  used  with  64  species  and  268  reactions  [19]. 
The  extended  Zeldovich  mechanism  was  used  for  NO  calculations.  In  order  to 
accommodate  the  entire  engine  load  spectrum,  the  drift  velocity  was  determined  using 
both  the  thermodynamic  conditions  and  the  electric  field  [19].  For  these  tests,  the  electric 
field  was  80  kV/m.  Slightly  different  expressions  result  for  the  ion  current  depending  on 
which  force  dominates  [19].  In  the  region  where  both  forces  are  important,  a  weighted 
linear  combination  is  applied  [19].  Both  positive  and  negative  ions  are  examined  in  the 
chemical  kinetic  model.  Negative  ions  are  also  governed  by  Saha’s  equation  (5)  through 
the  example  reaction: 


M' +  E’ion  ^  M  +  e'  (11) 

where  the  negative  ion  (M')  is  in  a  ground  state  and  E’ion  is  the  energy  require  to 
neutralize  the  negative  ion  [19].  The  current  will  be  the  summation  of  electrons,  positive 
ions,  and  negative  ions  [19].  This  method  is  complicated  by  the  fact  that  the  ionization 
depends  upon  the  thermodynamic  state  that  the  sensor  is  trying  to  detect  [19].  Therefore 
additional  information  and  assumptions  are  required  to  solve  the  problem  [19].  Einally,  a 
range  of  engine  velocities,  throttle  positions,  torques,  ignition  timings,  and  lambda  are 
investigated  to  experimentally  validate  the  improved  model  [19]. 

Using  the  previous  model,  the  peak  values  for  ion  current  and  pressure  are 
correlated  and  the  correlation  coefficient  is  found  to  be  0.6  [19].  The  “not  particularly 
high”  [19]  coefficient  is  the  result  of  erratic  ion  currents  for  each  cycle.  Saitzkoff  et  al. 
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[19]  found  that  the  ion  current  decreased  at  lower  equivalence  ratios  when  results  were 
averaged  over  500  cycles.  At  leaner  mixtures  the  NO  concentration  increases  and 
therefore  the  ion  current  is  expected  to  increase  as  well.  The  reason  for  the  decreasing  ion 
current  at  leaner  mixtures  is  the  decrease  in  temperature  and  therefore  lower  ionization 
ratios  [19].  The  dominant  electronegative  species  was  found  to  be  the  hydroxyl  radical 
(OH)  [19].  Results  show  long  lived  hydrocarbon  species  are  not  important  to  the  ion 
current  [19].  Electrons  were  found  to  be  the  dominant  charge  carrier  due  to  the  higher 
drift  velocity  resulting  from  their  lower  mass  [19].  Saitzkoff  et  al.  also  found  a  strong 
correlation  in  time  between  the  maximum  peak  current  in  the  post-combustion  region  and 
the  maximum  peak  pressure  [19].  Slight  differences  were  explained  by  the  difference 
between  the  maximum  pressure  of  the  gas  and  the  maximum  density  [19].  Over  a  large 
number  of  cases  the  correlation  coefficient  of  0.8  was  obtained  without  filtering  of  the 
current  [19].  Although  low  load  driving  conditions  had  low  correlations,  in  averaged 
cases  the  correlation  of  predicted  pressure  was  above  0.95  [19].  Overall  Saitzkoff  et  al. 
showed  that  pressure  can  be  predicted  by  the  ion  current  although  not  currently  as 
accurate  as  desired. 

11,5  Shock  and  Detonation  Waves 

The  flame  structure  in  both  space  and  time  can  be  highly  turbulent  especially 
when  interacting  with  waves  [14].  This  complexity  was  examined  for  detonation  waves 
in  hydrogen-oxygen  mixtures  [14].  Results  for  detonation  waves  show  the  shock  has  little 
influence  on  the  ionization  compared  to  the  flame  or  the  detonation  wave  [14].  In  a 
hydrogen-oxygen  mixture  with  1.0  percent  N2  the  conductivity  was  found  to  be  4.5  x  10'"^ 
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and  4  x  10'^  (ohm  cm)'^  for  the  detonation  and  flames  respectively  [14],  Conversely,  the 
shock  conductivity  was  several  orders  of  magnitude  lower  at  5.4  x  10'  (ohm  cm)'  [14], 
Combustion  will  therefore  produce  ionization  levels  nine  orders  of  magnitude  larger  than 
a  shock  alone.  The  shock  ionization  was  attributed  to  purely  thermal  equilibrium  values 
[14].  NO  also  dominates  the  ion-producing  species  behind  detonations  [14].  Some 
impurities,  however,  could  become  highly  ionized  and  obscure  the  ionization  distinctions 
between  flames,  detonations  and  shocks  [14].  Assuming  effects  of  impurities  are 
insignificant,  the  generation  of  ions  by  the  shock  can  be  considered  negligible.  The 
ionization  behind  the  detonation  wave  should  also  be  larger  than  ordinary  flames. 

II.6  Structure  of  the  Pulse  Detonation  Wave 

While  the  physical  time  and  space  structure  of  a  detonation  wave  is  extremely 
complex  three  dimensional  phenomena,  the  one  dimensional  theoretical  structure  as 
described  by  Zeldovich,  von  Neumann,  and  Doring,  referred  to  as  ZND  wave  structure 
[6],  is  relatively  straightforward  as  shown  in  Figure  1. 
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Von  Neumann 
Pressure  Spike 


Figure  1.  ZND  Detonation  wave  structure  moving  from  left  to  right 


The  detonation  wave  begins  near  the  closed  end  of  the  tube  and  travels  through  the 
unbumed  fuel-air  mixture  towards  the  open  end  of  the  tube.  The  fuel-air  mixture  first 
encounters  a  strong  shock  wave  that  compresses  the  mixture,  elevating  both  the 
temperature  and  pressure  [6].  After  a  short  ignition  delay,  or  induction  zone,  the 
combustion  initiates  creating  Rayleigh  type  heat  addition  into  the  flow  [13].  Once  the 
combustion  is  complete  at  state  2,  the  flow  is  at  Chapman-Jouguet  conditions  for  a  self- 
sustaining  detonation  [6].  The  shock  wave  and  combustion  zone  are  closely  coupled 
phenomena  necessary  for  a  detonation  wave.  The  speed  of  the  wave  relative  to  the  burned 
mixture  is  sonic,  whereas  relative  to  the  unburned  mixture  the  wave  is  supersonic.  The 
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speed  of  the  wave  is  driven  by  the  heat  release  rate  of  the  eombustion.  A  large  pressure 
spike,  termed  the  von  Neuman  spike,  oecurs  immediately  behind  the  shoek  wave  [6].  As 
the  wave  passes,  the  flow  expands  to  the  level  of  the  elosed  end  through  rarefaetion 
waves  [6]. 

The  total  width  of  the  shoek,  induetion  zone,  and  eombustion  zone  is  on  the  order 
of  one  eentimeter  [6].  Measured  wave  speeds  of  1957  meters  per  seeond  (m/s)  were 
produeed  from  the  Air  Foree  Researeh  Lab,  Propulsion  Direetorate  PDF  (AFRL/PRTS) 
[20].  Therefore,  the  time  for  the  detonation  wave  to  travel  a  1  meter  tube  is  on  the  order 
of  0.5  ms.  The  time  for  the  entire  detonation  wave  strueture,  assuming  a  length  of  1  em, 
to  pass  a  point  in  the  tube  is  on  the  order  of  15  mieroseeonds  (ps)  or  66  kilohertz  (kHz) 

o 

[6].  The  shoek  wave  itself  has  a  length  of  several  moleeular  mean  free  paths  (6  X  10'  m) 
resulting  in  a  time  on  the  order  of  30  pieoseeonds  (fs)  or  33  gigahertz  (GHz)  to  pass  a 
point  in  the  tube.  Assuming  the  time  for  expansion  is  larger  than  the  detonation  wave 
travel  time,  the  time  for  the  pressure  to  deeay  from  state  2  to  the  elosed  end  wall  is  on  the 
order  of  1  ms. 

IL7  Comparison  of  Engine  Times 

IC  engines  and  the  PDF  operate  on  different  time  seales.  Depending  on  the  speed 
of  the  IC  engines,  the  eombustion  and  post-eombustions  time  is  on  the  order  of  20  ms 
whereas  the  PDF  is  on  the  order  of  1  ms.  The  PDF  is  roughly  one  order  of  magnitude 
faster  than  IC  engines.  Laminar  flame  speed  in  air,  for  eomparison,  is  on  the  order  of  0.7 
m/s  for  hydroearbons  and  0.4  m/s  for  methane  [12:130-131].  Detonation  waves  are 
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therefore  three  and  a  half  orders  of  magnitude  faster  than  laminar  flame  speeds.  These 
different  time  scales  play  an  important  role  in  ionization  and  measurement  techniques. 

11,8  Expected  Ion  Current  in  the  PDE 

The  initial  ionization  level  is  extremely  small  in  low  temperature  regions  such  as 
the  unburned  mixture  in  the  PDE  before  the  arrival  of  the  detonation  wave.  The  shock 
wave  increases  the  degree  of  ionization  through  thermal  heating  and  remains  roughly 
constant  through  the  induction  zone.  As  previously  discussed,  the  combustion  zone 
sharply  increases  the  ionization  level  by  nine  orders  of  magnitude  higher  than  the  shock 
wave,  due  to  both  chemi-ionization  and  thermal  effects  at  elevated  temperatures.  This 
high  super-equilibrium  ionization  level  will  then  decay  down  to  equilibrium  levels  at 
lower  temperature. 

The  low  ionization  levels  in  the  unbumed  fuel-air  mixture  will  make  the 
measurement  of  an  ion  current  above  the  noise  extremely  difficult.  The  first  increase  in 
ion  current  would  be  created  from  the  shock.  The  short  lengths  of  the  shock  and  induction 
periods  make  measuring  this  current  impractical.  The  combustion  zone  will  create  a  large 
increase  in  the  ion  current.  The  ion  current  decay  rate  from  this  super-equilibrium  level  is 
a  function  of  the  square  of  pressure  through  molecular  recombination.  The  pressure 
during  the  expansion  region  could  be  determined  by  measuring  the  ion  decay  rate  and 
correlating  it  to  the  pressure. 
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11.9  Determining  Pressure  from  the  Ion  Current  in  the  PDE 

This  method  could  provide  useful  pressure  information  between  state  2,  as  shown 
in  Figure  1,  and  the  closed  end  wall.  The  von  Neumann  pressure  spike  would  be 
impractical  to  measure  due  to  the  short  duration  and  low  ionization  level  compared  to  the 
end  of  the  combustion  zone.  Because  this  method  relies  upon  the  ion  decay  and  not 
equilibrium  concentrations,  the  method  can  only  be  applied  to  unsteady  processes  and  not 
steady  state  conditions  typically  measured  by  conventional  piezoelectric  pressure 
transducers.  Like  the  IC  engines,  the  ion  current  in  the  PDE  can  also  be  measured  using  a 
spark  plug.  The  durability  of  the  spark  plug  is  inherently  suited  for  the  harsh  environment 
of  the  PDE. 

The  ion  current  in  the  PDE,  unlike  the  IC  engines,  should  not  experience  a  second 
ionization  peak  due  to  NO  formation.  Although  NO  is  still  the  dominate  ion  producing 
species  behind  detonation  waves,  the  relative  concentration  is  lower  than  in  IC  engines. 
The  order  of  magnitude  quicker  processes  of  the  detonation  wave  and  expansion  zone 
provide  less  time  for  the  slow  chemical  kinetics  to  form  NO.  Eurther,  the  lower  pressures 
in  the  PDE  will  lower  the  partial  pressure  and  thus  density  of  any  NO  produced.  Chemi- 
ionization  will  be  the  dominant  factor  in  ion  production  for  the  PDE. 

11.10  Derivation  of  the  Pressure  and  the  Ion  Current  Decay  Rate  Relationship 

A  simple  derivation  can  show  how  the  decay  rate  of  the  ion  current  relates  to  the 
pressure.  The  electric  current  /,  by  definition,  is  the  rate  that  charge  passes  through  a 
surface.  In  this  case  the  surface  is  the  exposed  area  of  the  ion  probe.  Current  can  also  be 
expressed  in  the  form: 
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/  =  nAv^  q 


(12) 


where  n  is  the  number  density  of  the  eharge  carriers,  A  is  the  surface  area  of  the  probe, 
is  the  drift  velocity  and  q  is  the  charge.  Assuming  the  net  motion  is  due  to  the  electric 
field,  an  ion  can  be  accelerated  by  the  electric  field  until  it  collides  with  another 
molecule.  The  average  speed  of  the  electrons  is  not  considered  because  the  velocity  after 
collision  is  randomly  directed  and  does  not  contribute  to  the  drift  velocity.  The  drift 
velocity  is  then  expressed  by: 


q  Et 
m 


(13) 


where  E  is  the  electric  field,  m  is  the  mass  of  the  electron,  and  x  is  the  mean  free  time. 
The  mean  free  time  is  simply: 


r  =  ■ 


(14) 


where  ^  is  the  mean  free  path  and  Vavg  is  the  average  electron  velocity.  The  drift  velocity 
is  therefore: 
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The  mean  free  path  is  proportional  to  temperature  divided  by  pressure  and  the  average 
velocity  is  proportional  to  temperature 


Zoz-  (16) 

P 
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where  T  is  the  temperature  and  P  is  the  pressure.  The  drift  velocity  is  therefore 
proportional  to  the  electric  field  divided  by  the  pressure. 
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The  drift  velocity  is  also  a  weak  function  of  temperature.  Since  it  is  not  a  dominate  term, 
it  may  be  neglected  for  these  purposes.  The  change  in  current  with  respect  to  time  is: 


dt 


(19) 
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As  previously  discussed,  the  ion  decay  rate  is  a  function  of  the  square  of  pressure  due  to 
the  dominant  bi-molecular  recombination.  Assuming  a  constant  electric  field,  the  ion 
current  decay  rate  is  proportional  to  the  pressure: 


dl 

—  oc 
dt 


rn 
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(20) 


From  this  relationship,  the  pressure  can  be  determined  by  measuring  the  ion 
current  decay  rate.  By  using  the  spark  plug  to  measure  the  ion  current  and  applying  this 
simple  relationship,  a  complementary  pressure  sensing  technique  can  be  developed. 
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Ill,  Experimental  Approach 


111.1  Chapter  Overview 

This  chapter  describes  how  the  predieted  behavior  of  the  previous  ehapter  will  be 
tested.  The  test  setup  is  explained  along  with  instrumentation  and  data  flow. 

111.2  Combustion  Bomb 

Ultimately  the  spark  plug  is  sought  to  be  an  additional  pressure  sensing  deviee  in 
the  PDE.  The  unsteady  nature  of  the  PDE  and  the  harsh  environment  make  eorrelations 
between  pressure  and  ion  eurrent  diffieult.  Sinee  the  PDE  elosely  follows  a  eonstant 
volume  eyele,  a  eombustion  bomb  experiment  ean  allow  investigation  into  the  ion  eurrent 
dependenee  upon  pressure.  Conventional  piezoeleetrie  pressure  transducers  ean  be  used 
in  a  eombustion  bomb  with  high  aecuraey.  The  eorrelations  developed  in  the  eonstant 
volume  proeess  can  be  applied  to  the  PDE  with  minimal  modifieations. 

Eventually  the  PDE  is  desired  to  run  on  praetieal  hydroearbon  fuels.  Eor  ease  of 
use,  methane  will  be  used  as  the  fuel  in  the  eombustion  bomb.  Dry  air  will  used  for  the 
oxidizer  for  both  ease  of  use  and  elose  approximation  to  PDE  operating  conditions. 

While  the  ehemistry  of  methane-air  reaetions  ean  vary  from  hydroearbon-air 
eombinations,  the  methane-air  mixture  allows  an  easy  first  investigation. 

The  laminar  methane-air  flame  speed  within  the  eombustion  bomb  is 
approximately  three  orders  of  magnitude  slower  than  the  detonation  wave  as  previously 
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discussed.  The  laminar  flame  speed  allows  sharper  time  resolution  into  the  proeess 
without  expensive  high  speed  instrumentation  needed  for  a  detonation  wave. 

This  experiment  is  designed  around  a  one  half  liter  stainless  steel  pressure  vessel 
rated  to  2000  pounds  per  square  ineh  (psi)  as  shown  in  Figure  2. 


Figure  2.  Overview  of  Experiment  Configuration 


A  stainless  steel  lid  with  six  aeeess  ports  seals  the  7.62  em  (3  in.)  outer  diameter  (OD) 
pressure  vessel.  Three  of  the  ports  are  1.429  eentimeters  (em)  or  9/16**^  ineh  (in.)  in 
diameter  while  the  other  three  symmetrie  ports  are  1.111  em  (7/1 6'*^  in.)  in  diameter.  A 
grounding  point  is  loeated  in  the  eenter  of  the  lid.  The  depth  of  all  of  the  ports  is  3.175 
em  (1.25  in.). 
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III,3  Ion  Sensor 


Three  different  ion  sensors  were  ereated  all  based  on  a  Champion  RC12LYC 
spark  plug  used  in  multiple  automotive  engines  as  shown  in  Figure  3.  The  plug  has  a 
measured  resistanee  of  approximately  57.5  kH.  This  spark  plug  was  seleeted  because  of 
the  low  cost  and  long  center  electrode.  The  Champion  spark  plug  also  had  a  compression 
washer  needed  for  a  tight  seal  against  the  un-tapered  top  of  the  vessel  lid. 


Figure  3.  Champion  RC12LYC  Spark  Plug 


111,2,1  Short  Probe 

The  first  ion  sensor,  referred  to  as  the  short  probe,  simply  had  the  side  prong 
removed  so  only  the  center  electrode  was  exposed  as  shown  in  Figure  4. 
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Figure  4.  Short  Probe  with  side  prong  of  RC12LYC  spark  plug  removed 

This  allowed  electrons  more  direct  access  to  the  center  electrode  thereby  increasing  the 
strength  of  the  ion  current.  The  0.257  cm  diameter  center  electrode  is  0.884  cm  long  with 
the  bottom  0.122  cm  exposed  and  the  upper  portion  covered  with  a  0.762  cm  diameter 
insulating  ceramic  material  originally  part  of  the  spark  plug.  When  placed  in  a  port  in  the 
lid  of  the  vessel,  the  bottom  tip  was  recessed  in  to  the  threaded  port  by  0.57  cm.  Centered 
in  the  port,  the  side  of  electrode  is  0.47  cm  inch  from  the  threaded  wall  of  the  port. 

111,3,2  Medium  Probe 

The  second  ion  sensor,  referred  to  as  the  medium  probe,  also  had  the  side  prong 
removed  but  the  center  electrode  was  extended  to  2.54  cm  (1.0  in.)  as  shown  in  Figure  5. 
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Figure  5.  Medium  Probe  with  2.54  cm  center  electrode 


This  extension  was  created  by  removing  some  of  the  original  ceramic  material  and 
connecting  the  original  center  electrode  to  a  2.54  cm  (1  in.)  steel  extension  using  a  0.635 
cm  (0.25  in.)  OD  steel  covering  with  a  small  set  screw.  The  electrode  extension  had  a 
diameter  of  0.267  cm,  similar  to  the  original  electrode  with  a  0.257  cm  diameter.  When 
placed  into  the  lid,  the  tip  of  the  medium  probe  extended  1.524  cm  below  the  bottom 
surface  of  the  lid  and  1 .27  cm  from  the  vessel  side  wall. 

IIL4,3  Long  Probe 

The  third  ion  sensor,  referred  to  as  the  long  probe,  also  had  the  side  prong 
removed  but  the  center  electrode  was  extended  by  10.16  cm  (4.0  in.)  as  shown  in  Figure 

6. 
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Figure  6.  Long  Probe  with  10.16  cm  insulated  extension 


A  0.3175  cm  (0.125  in.)  inner  diameter,  0.635  cm  (0.25  in.)  OD  insulating  ceramic  tube 
was  placed  over  the  extending  electrode  leaving  0.635  cm  (0.25  in.)  of  the  probe  exposed. 
A  larger  1.08  cm  OD,  0.80  cm  ID  ceramic  tube  covered  the  connecting  section  of  the 
electrode.  Both  ceramic  tubes  have  a  high  electrical  resistance  and  were  secured  to  the 
electrode  using  blue  RTV  silicon  designed  for  automotive  applications.  Ultra  high 
temperature  RTV  was  not  used  because  of  the  conducting  properties  of  the  copper 
additive.  When  placed  into  the  vessel  lid  the  tip  of  the  electrode  extended  9.83  cm  below 
the  lid.  This  location  of  the  tip  also  corresponds  to  a  distance  of  approximately  5.0  cm 
from  the  bottom  of  the  vessel  and  approximately  1 .27  cm  from  the  vessel  side  wall.  The 
distance  from  the  side  of  the  vessel  may  vary  slightly  since  the  extension  was  not 
attached  perfectly  straight.  The  distance  from  the  tip  to  the  wall  could  vary  by  0.12  cm 
depending  on  the  final  rotation  of  the  probe  into  the  lid. 

III,4  Fuel  and  Air  System 

Both  the  methane  and  the  dry  air  entered  the  vessel  through  one  small  port  in  the 
lid.  This  port  also  acted  as  the  exit  port  for  the  combusted  products.  Plumbing  for  the  dry 
air  and  methane  was  accomplished  using  0.635  cm  (0.25  in.)  soft  copper  tubing  rated  to 
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225  psi.  Brass  Swagelok  provided  easy  and  reliable  eonneetions  for  the  eopper  tubing. 
The  brass  valves  were  rated  to  3000  psi.  A  bloek  diagram  of  the  air  and  methane  fuel 
system  is  shown  in  Figure  7. 


Figure  7.  Fuel  and  air  system  bloek  diagram 
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111,4.1  Methane 


The  source  of  methane  was  a  commercial  grade  standard  K-bottle  mounted  to  the 
lab  table.  The  pressure  of  the  methane  was  initially  controlled  by  a  regulator  on  the 
bottle.  The  regulator  was  then  connected  to  another  valve  using  copper  tubing.  A  check 
valve  then  ensured  proper  flow  of  the  methane  and  eliminated  any  backflow  concerns.  A 
needle  valve  then  tightly  controlled  the  pressure  of  the  methane.  Before  the  methane 
entered  the  vessel,  another  shut-off  valve  allowed  the  mixture  within  the  combustion 
bomb  to  be  closed  off  from  the  other  plumbing. 

IIL4,2  Dry  Air 

Two  separate  sources  of  dry  air  were  available  for  the  experiment  also  shown  in 
Figure  7.  A  Jun-Air  model  3-1.5  air  compressor  provided  a  120  psi  source  of  air.  A  water 
and  oil  separator  attached  to  the  compressor  ensured  low  levels  of  humidity  and 
contaminants  although  exact  levels  were  not  measured.  A  second  600  psi  source  of  dry 
air  was  brought  into  the  lab  from  an  outside  tank.  This  second  source  allowed  for  test 
cases  above  8  atm  although  not  required.  Both  sources  of  dry  air  connected  to  valves  and 
then  to  a  T  connector.  Similar  to  the  methane  source,  a  check  valve  and  needle  valve 
properly  controlled  the  direction  and  pressure  of  the  dry  air  entering  the  pressure  vessel. 

IIL4,3  Vacuum  Pump 

A  Franklin  Electric  0.5  horsepower  vacuum  pump  was  connected  to  the  vessel 
after  the  control  valve  with  plastic  tubing  since  high  pressures  would  not  be  seen  by  the 
pump. 
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111.4,4  Exhaust 


The  combustion  products  were  removed  from  the  vessel  through  plastic  tubing 
leading  to  an  exhaust  fan,  dumping  the  products  outside  the  building. 

111,4,5  Partial  Pressure  Control 

The  method  of  partial  pressures  ensured  accurate  control  of  the  equivalence  ratio. 
When  filling  the  vessel  with  methane  and  dry  air  the  vessel  pressure  was  measured  by  an 
Endevco  15  psi  absolute  (psia)  conventional  pressure  transducer.  The  calibrated  accuracy 
of  the  sensor  was  1%.  A  separate  Endevco  4428A  conditioning  box  powered  the  sensor 
and  controlled  the  calibration.  The  sensor  was  connected  to  the  vessel  through  a  small 
port  that  was  split  by  a  T  connector.  On  one  side  of  the  “T”  was  the  main  pressure 
transducer  used  during  the  experiment  to  correlate  the  ion  current.  On  the  other  side  of 
the  “T”  was  the  Endevco  transducer  separated  by  a  shut-off  valve.  When  filling  the 
vessel,  this  valve  was  open  to  allow  the  Endevco  sensor  to  accurately  measure  the  partial 
pressures.  During  the  experiment  and  any  other  times  where  pressures  were  above  15 
psia,  this  valve  was  closed  to  prevent  any  damage  to  the  sensor. 

The  Endevco  pressure  transducer  also  aided  with  the  vacuum  pump.  The  sensor 
ensured  consist  vacuum  levels  prior  to  introducing  the  methane  and  dry  air.  The  exact 
vacuum  level  the  pump  was  capable  of  was  unknown,  but  the  accuracy  of  the  Endevco 
sensor  ensured  a  vacuum  level  of  0.15  psia  or  lower. 
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111.4,6  Fuel  &  Air  Procedure 


The  vessel  was  first  pumped  down  to  an  assumed  level  of  0.15  psia.  Methane 
introdueed  to  the  system  inereased  the  pressure  until  the  desired  methane  partial  pressure 
was  aehieved.  Onee  the  methane  flow  stopped,  dry  air  was  added  to  the  vessel  to  the 
desired  total  pressure.  If  the  desired  total  pressure  was  14.7  psia  or  1  atmosphere  (atm), 
the  Endeveo  sensor  was  used  to  measure  the  vaeuum  level,  pressure  while  adding 
methane,  and  the  pressure  while  adding  dry  air.  If  the  desired  total  pressure  was  above  1 
atm,  the  valve  before  the  Endeveo  sensor  was  elosed  after  the  methane  was  added.  The 
dry  air  was  then  added  and  eontrolled  using  the  main  experiment  pressure  transdueer. 
After  the  individual  test  eompleted,  the  produets  were  sueked  out  the  exhaust  fan.  The 
system  was  then  flushed  with  dry  air  three  times  to  help  remove  any  eontaminants.  The 
vaeuum  pump  then  removed  any  remaining  eontaminants  and  proeedure  was  repeated  for 
additional  tests. 

III,5  Ignition  System 

The  fuel  and  air  mixture  is  ignited  using  a  traditional  automotive  induetive 
diseharge.  An  unmodified  Champion  RC12EYC  spark  plug,  as  shown  in  figure  3, 
produees  the  ignition  spark.  A  MSD  Blaster  3  (MSD-8223)  ignition  eoil  generates  a 
maximum  45,000  volts  to  the  spark  plug.  This  eoil  uses  a  tall  tower  to  improve  the  spark 
isolation  and  eoil  wire  attaehment.  The  reeommended  0.8  Q.  ballast  resistor  was 
eonneeted  in  between  the  ignition  eoil  and  the  power  supply.  An  HP  6033A  power  supply 
ereates  a  elean  high  eurrent  12  V  DC  supply  for  the  ignition  eoil.  The  6033A  has  a  RMS 
noise  level  of  3  milli-Volts  (mV).  Typieal  12  V  automotive  power  supplies  and  battery 
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chargers  easily  provide  the  required  power  for  the  ignition  but  also  generate  unacceptable 
noise  in  the  system.  This  occurs  because  the  ignition  spark  plug  and  the  ion  probe  share 
the  common  and  ground  signals.  The  ignition  coil  amplifies  any  noise  in  the  12  V  supply. 
A  low  noise  power  supply  such  as  the  HP  6033A  is  critical  to  reducing  overall  system 
noise. 

The  electrical  ignition  circuit  is  the  same  as  in  older  automotive  applications  as 
shown  in  Figure  8. 


Figure  8.  Ignition  Circuit 
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The  12  V  DC  flows  through  the  ballast  resistor  and  then  the  eoil,  generating  a  strong 
magnetie  field.  The  eurrent  then  flows  through  a  meohanieal  40  Amp  rated  relay  to  the 
power  supply  eommon  and  ground.  The  position  of  the  relay  was  eontrolled  by  a 
normally  elosed  single  pole  double  throw  momentary  switeh.  The  12  V  DC  powered 
relay  was  normally  elosed  until  the  momentary  switeh  was  pressed  eausing  the  relay  to 
break  the  ignition  cireuit  foreing  the  high  voltage  diseharge.  The  relay  was  better  suited 
than  the  switeh  to  break  the  eireuit  beeause  of  the  quiek  break  of  the  eonnection  and 
higher  durability.  The  swift  break  of  the  eireuit  is  erueial  to  generating  a  strong  spark 
aeross  the  spark  plug.  When  the  eireuit  is  broken  by  the  relay,  the  strong  magnetie  field 
induees  a  high  voltage  diseharge  through  the  seeondary  windings  of  the  eoil.  An 
automotive  0.26  miero-farad  eapaeitor  was  eonneeted  in  parallel  to  the  relay  and  ground. 
This  eapaeitor  forees  the  potential  at  both  sides  of  the  relay  to  remain  at  ground  potential 
also  eontributing  to  a  elean  break  of  the  eireuit.  Without  the  eapaeitor  at  the  relay,  the  eoil 
will  not  produee  the  required  high  voltage.  The  high  voltage  aeross  the  spark  plug  eauses 
the  mixture  to  breakdown  into  plasma  ereating  a  short  but  high  energy  region  that  forees 
ignition  in  the  rest  of  the  fuel-air  mixture. 

TTT.6  Instrumentation 

In  order  to  investigate  the  relationship  between  the  ion  eurrent  and  the  pressure 
several  sensors  are  required. 


37 


IIL6.1  DC  Voltage  for  Ion  Probes 

The  ion  probes,  as  previously  deseribed,  eapture  the  ionization  levels  through  a 
eurrent  aeross  the  probe  when  energized  by  a  DC  voltage.  Eaeh  probe  was  energized  into 
an  ion  sensor  by  applying  a  positive  lOV  DC  from  a  ME  83B829  power  supply.  The 
voltage  was  held  within  a  toleranee  of  0.01  V.  The  positive  voltage  was  plaeed  on  the 
eenter  eleetrode  while  the  pressure  vessel  was  the  eommon  side  of  the  signal.  Other  DC 
voltages  of  positive  5,  20,  and  40  volts  were  also  investigated. 

111,6,2  Current  Measurement 

A  Keithley  6487  pieoammeter  measures  the  eurrent  aeross  the  probe.  This 
pieoammeter  has  a  measuring  rate  of  1000  Hz.  Eor  the  experiments,  neither  the 
dampening  funetion  nor  the  internal  voltage  souree  were  used.  Two  different  seales  were 
manually  set  for  testing:  20  mieroamps  (pA)  and  2  pA.  The  internal  buffer  and  ability  to 
eommand  the  unit  were  also  not  used.  Instead,  a  high  speed  data  aequisition  eard  (DAC) 
direetly  eaptured  the  analog  output  of  the  meter. 

Shown  in  Eigure  9,  the  meter  was  plaeed  in  the  ion  probe  eireuit  between  the 
power  supply  and  the  probe  sinee  the  meter  must  be  in  series  with  the  desired  eurrent. 
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Figure  9.  Ion  Probe  Cireuit 


This  picoammeter  also  required  the  input  high  end  of  the  device  to  be  connected  to  the 
high  resistance  side  of  the  circuit.  In  this  case  the  57.5  kH  resistance  built  into  the 
Champion  RC12LYC  became  the  high  resistance  portion  of  the  circuit.  Since  the  polarity 
of  the  power  supply  and  the  picoammeter  are  opposite,  the  current  measured  across  the 
picoammeter  will  be  displayed  as  negative  although  physically  the  opposite  is  true. 

To  summarize  the  ion  probe  circuit,  the  ME  83B829  power  supply  generates  a 
small  constant  voltage  and  corresponding  current.  This  current  flows  through  the 
picoammeter  and  then  flows  across  the  ion  probe.  The  common  side  of  the  probe  is  the 
thread  of  the  spark  plug  that  is  mounted  into  the  vessel.  The  ground  point  of  the  vessel  is 
connected  to  a  system  ground  point  that  is  connected  to  the  return  side  of  the  power 
supply. 
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111.6,3  Pressure  Sensors 


For  pressure  measurements  two  different  Omega  transdueers  were  used 
depending  on  the  maximum  expeeted  pressure  in  the  vessel.  The  eonventional  Omega 
PX303-300A5V  transdueer  measured  pressures  for  test  eases  with  an  initial  pressure  of  3 
atm  or  less.  This  transdueer  useful  range  is  from  0  to  300  psia  with  a  response  time  of  1 
ms.  For  test  eases  above  an  initial  pressure  of  3  atm,  a  similar  Omega  PX303-1KG5V 
transdueer  measured  pressures  from  0  to  1000  psi  gauge  (psig).  The  gauge  pressure 
readings  were  eonverted  into  absolute  pressures  by  adding  the  atmospherie  pressure  in 
the  lab  as  measured  by  a  Druek  DPI- 141  digital  barometer.  While  the  aeeuraey  of  both 
transdueers  is  the  same  pereentage,  the  higher  range  of  the  1000  psig  transdueer  provides 
less  resolution  into  the  pressure.  An  Omega  PSS-15  power  supply  powered  the  Omega 
transdueers.  Both  transdueers  produeed  a  0.5-5. 5  V  signal  eorresponding  to  the 
minimum  and  maximum  pressure  respeetively. 

111,6,4  Thermocouples 

Two  K-type  thermocouples  were  used  in  the  experiment  to  measure  the  ambient 
temperature  in  the  lab  and  the  internal  temperature  of  the  combustion  bomb.  The  ambient 
sensor  was  directly  connected  to  a  SCXI-1 1 12  signal  conditioning  module  described  later 
in  this  chapter.  The  internal  thermocouple  was  inserted  into  the  vessel  through  one  of  the 
small  ports  in  the  lid  and  secured  with  a  graphite  compression  fitting.  The  tip  of  the 
thermocouple  extended  approximately  0.32  cm  below  the  bottom  surface  of  the  lid.  This 
position  was  selected  to  measure  gas  temperature  without  interfering  with  other  sensors. 
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111.6,5  Band  Heater 


A  Chromalox  band  heater,  loeated  at  the  bottom  the  pressure  vessel,  inereases  the 
temperature  of  the  gases  in  the  lower  portion  of  the  vessel.  Due  to  the  generated 
buoyaney,  the  dense  gases  rise  towards  the  top  of  the  vessel  where  the  eooler 
temperatures  allow  density  to  inerease  eausing  the  gas  to  drop  baek  towards  the  bottom. 
This  buoyaney  inereases  the  mixing  of  the  lighter  methane  with  the  denser  dry  air.  In 
addition  to  the  buoyaney  effeet,  the  band  heater  also  raises  the  overall  temperature  of  the 
pressure  vessel.  The  band  heater  is  eontrolled  by  an  Omega  CSC32  beneh  top  eontroller 
that  allows  the  temperature  to  be  manually  set  or  remotely  through  a  computer  program 
using  the  serial  port.  The  controller  measures  the  temperature  of  a  thermocouple  inserted 
in  between  the  pressure  vessel  and  the  band  heater.  This  indicted  the  temperature  of  the 
band  heater  and  not  the  internal  temperature  of  the  gases  in  the  vessel  nor  the  ambient 
temperature. 

The  120  V  AC  powered  resistance  band  heater  operated  intermittently.  A 
proportional-integral-derivative  (PID)  controller  maintained  the  set  temperature  for  the 
heater.  The  controller  was  tuned  using  the  software  included  with  the  controller.  Using 
this  controller,  the  temperature  of  the  band  heater  stayed  within  2  °C  of  the  set  point  at  all 
times.  In  addition  to  displaying  the  current  temperature,  the  PID  controller  also  indicted 
when  current  flowed  through  the  resistance  heater. 

III,7  Data  Flow 

The  data  from  the  sensors  traveled  through  one  of  two  routes  as  indicted  in  Figure 

10. 
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Figure  10.  Experiment  Data  Flow 


The  pieoammeter  measured  the  ion  eurrent  and  sent  the  results  to  a  high  speed  DAC 
where  a  Lab  View  program  eaptured  the  data.  The  pressure  sensors  and  thermoeouples 
used  a  SCXl-1000  ehassis  to  pass  the  data  before  being  read  by  a  slower  DAC  again 
where  a  LabView  program  eaptured  the  data.  The  data  was  then  plotted  and  analyzed 
using  MatLab. 

IIL7.1  SCXI-1000  Chassis 

The  National  Instruments  (NI)  SCXI-1000  ehassis  provides  an  easily  eonfigurable 
platform  to  pass,  eondition,  and  multiplex  input  and  output  signals.  The  ehassis  has  spaee 
for  up  to  four  modules  although  only  three  were  used  for  this  experiment.  A  SCXI-1 112 
module  was  used  to  amplify  and  eondition  both  the  internal  and  ambient  thermoeouple 
signals.  A  SCXI-1 180  module  provided  a  direet  pass  through  of  signals  straight  to  the 
DAC.  A  SCXI-1302  breakout  box  eonneeted  to  the  SCXI-1 180  module  allowed  simple 
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hookups  to  the  50  pin  feed-through  panel.  The  signals  from  the  Omega  pressure 
transducers  were  connected  to  the  SCXI-1302  breakout  box  since  no  conditioning  of  the 
signal  was  required.  The  previously  described  double  throw  ignition  switch  used  one 
throw  to  break  the  ignition  circuit  causing  the  high  voltage  discharge.  The  other  throw  on 
the  switch  changed  a  normally  0  V  signal  to  a  +5  V  DC  signal  that  was  read  by  both  the 
slow  and  fast  DACs.  This  timing  signal  triggered  the  data  capture  of  the  combustion 
event.  A  third  SCXI-1 124  digital  to  analog  converter  module  with  a  SCXI-1325  breakout 
box  provided  analog  out  signals  to  control  other  hardware.  The  speed  of  SCXI-1 000 
chassis  has  a  maximum  200  kilo-Samples  per  second  (kS/s)  for  all  channels.  The  signals 
routed  through  the  SCXI-1302  module  were  not  part  of  the  SCXI-1000  sample  rate  but 
rather  the  rate  of  the  DAC  card  itself  In  this  configuration  the  SCXI-1000  required  an  E- 
series  NI  DAC. 

IIL7.2  NI6024EDAC 

A  NI  6024E  DAC  provided  data  capture  for  up  to  eight  12-bit  channels  at  200 
kS/s.  This  NI  E-series  card  was  able  to  use  the  updated  EabView  mx  drivers  in  addition 
to  the  traditional  EabView  drivers.  This  card,  however,  does  not  support  hardware 
triggering.  This  slower  DAC  card  captured  the  pressure  and  thermocouple  data.  The  +5  V 
ignition  timing  signal  was  also  captured  by  this  card.  The  card  also  read  the  multiplexed 
signal  from  the  SCXI-1000  chassis. 


43 


IIL7.3  NI6110DAC 


A  NI  61 10  DAC  supplied  data  capture  on  four  12-bit  channels  of  up  to  5  mega¬ 
samples  per  second  (MS/s)  per  channel.  While  this  fast  data  capture  rate  was  excessive 
for  this  experiment,  the  extra  speed  could  be  useful  for  testing  in  the  PDE  itself.  This  card 
was  not  supported  by  the  new  LabView  mx  drivers  but  was  capable  of  hardware 
triggering.  The  ion  current  was  captured  by  one  channel  of  this  device  while  another 
channel  was  used  to  capture  the  +5  V  ignition  timing  signal.  This  card  does  not  support  a 
SCXl-1000  chassis. 

111.7,4  Labview  Program 

The  Labview  7.0  software  package  was  used  to  control  experiment  conditions, 
capture  sensor  data,  and  display  results.  LabView  operated  under  Windows  2000 
Professional  on  a  computer  with  a  1000  MHz  Pentium  processor  and  256  MB  of 
memory. 

111,7,4,1  Program  Modes 

A  Labview  Virtual  Instrument  (VI),  or  program  code,  was  written  to  operate  in 
two  different  modes.  The  first  mode  continuously  displayed  sensor  data  onto  the 
computer  screen  allowing  quick  checks  of  system  operation  and  proper  initial  conditions 
before  initiating  the  experiment.  Once  all  conditions  were  satisfactory,  the  program  was 
manually  changed  to  the  data  acquisition  mode.  Under  this  mode  the  LabView  program 
began  searching  the  incoming  data  for  the  ignition  timing  signal.  Once  the  timing  signal 
changed  from  0  to  +5  V,  the  data  acquisition  began.  Data  0.1s  before  and  1  s  after  this 
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trigger  were  transferred  from  the  DACs  internal  buffer  and  into  the  main  system  memory. 
Any  data  in  the  buffers  before  or  after  this  period  was  discarded.  With  the  useful 
information  in  memory,  the  data  was  formatted  and  saved  into  two  separate  data  files, 
one  for  each  DAC. 

111.7.4.2  DAC  Drivers 

Since  only  the  E-series  card  was  supported  by  the  new  LabView  mx  drivers,  the 
main  VI  had  to  be  written  using  the  older  traditional  Vis  and  routines.  This  complicated 
the  program  since  easier  improved  tasks  in  the  new  mx  driver  could  not  be  used.  Despite 
this  obstacle,  the  main  VI  worked  well  once  bugs  in  the  program  were  worked  out. 

111.7.4.3  Triggers 

The  lack  of  a  hardware  trigger  on  the  E-series  card  also  complicated  the 
programming.  The  fast  5  MS/s/channel  DAC  provided  hardware  trigger  capability  where 
the  card  itself  searched  the  incoming  data  for  the  required  trigger  value.  In  this  way, 
hardware  triggering  uses  minimal  computer  system  resources  and  is  time  accurate  within 
half  the  sample  time.  To  overcome  this  deficiency  with  the  E-series  card,  a  software 
conditional  trigger  scheme  was  applied.  This  scheme  transfers  the  data  from  the  internal 
buffer  on  the  slow  DAC  into  system  memory  where  the  EabView  program  searches  the 
data  for  the  desired  trigger  value.  This  conditional  software  trigger  can  take  considerable 
processing  power  and  system  memory  depending  of  the  data  rate.  Since  only  the  slow 
DAC  used  this  scheme,  the  burden  of  the  computer  did  not  reduce  the  functionality  of  the 
rest  of  the  EabView  program.  In  other  words,  the  computer  was  still  powerful  enough  to 
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complete  the  other  required  LabView  tasks  at  the  same  time  without  impacting  results. 
Another  disadvantage  to  the  conditional  trigger  was  the  lower  time  accuracy.  When 
transferring  the  data  from  the  internal  buffer  to  the  system  memory  there  is  a  greater 
ehanee  for  skewing  of  the  measurement  time.  The  times  between  samples  remains 
eonstant  but  some  error  could  be  introduced  that  varies  the  triggering  time.  Both  triggers, 
despite  reading  the  same  timing  signal  at  the  same  trigger  level,  could  have  slightly 
different  start  times.  Although  not  desirable,  this  timing  error  was  acceptable  for  this 
experiment. 

III.7.5  Acquisition  Rates 

Aequisition  rates  of  both  DACs  were  selected  to  make  the  instruments  the 
limiting  faetor  in  the  data  rate  and  not  the  cards  themselves.  A  rate  of  10  kS/s/channel 
was  seleeted  for  the  slow  E-series  DAC.  This  rate  was  easily  fast  enough  to  eapture  the 
pressure  and  temperature  and  provided  good  time  resolution  into  the  timing  signal 
without  overburdening  the  computer.  A  rate  of  100  kS/s/channel  was  selected  for  the  fast 
6110  eard.  This  speed  easily  eaptured  the  ion  current  signal  from  the  pieoammeter 
without  creating  exeessive  amounts  of  data.  As  previously  described,  the  data  was 
acquired  for  a  total  of  1.1  s  with  0.1  s  of  data  before  the  trigger. 

III.8  Shielding 

Initial  tests  eonducted  to  work  out  the  bugs  in  the  system  showed  high  levels  of 
noise  in  the  ion  eurrent.  Some  noise  was  expeeted  due  to  the  low  miero-amp  currents 
being  measured.  The  ion  eurrent  signal  was  analyzed  in  LabView  by  using  a  built-in  fast 
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Fourier  transform  (FFT)  and  then  plotted  the  results  against  the  frequeney.  This  revealed 
a  spike  in  the  signal  at  60  Hz.  This  noise  was  assumed  to  be  generated  from  a 
combination  of  a  poor  facility  ground,  fluorescent  lighting,  various  power  supplies, 
measurement  instrument  themselves,  and  the  nearby  computer.  Moving  the  setup  to  a 
lower  noise  electronics  lab  was  not  possible  due  the  hazardous  methane.  Instead,  the 
noise  was  reduced  by  creating  a  better  single  point  system  ground  and  improving 
shielding  on  signal  wires.  The  signal  wires  already  had  standard  wire  insulation  and  in 
some  cases  such  as  the  picoammeter  cable,  a  shield  protected  portions  of  the  wiring. 
Aluminum  tape  was  wrapped  around  the  unshielded  wires  and  then  connected  to  the 
appropriate  ground  point.  By  adding  the  shielding,  the  60  Hz  was  reduced  by  over  10  dB 
but  still  remained  the  dominate  source  of  noise  as  expected.  While  noise  is  not  desirable, 
the  ion  sensor  will  also  experience  noise  when  employed  in  the  PDF.  The  more  noise  the 
sensor  is  able  to  handle  in  the  laboratory  setting  the  easier  it  will  be  to  successful  employ 
the  sensor  in  the  PDE. 

III,9  Steady  State  Ion  Current  Experiment 

Initially,  the  steady  state  ion  currents  were  desired  to  be  a  useful  comparison 
against  levels  measured  during  and  after  combustion.  This  first  experiment  used  a 
Classman  model  PS/FC20R06.0-1 1  high  voltage  power  supply  capable  of  up  to  20,000 
volts  positive  or  negative.  This  power  supply  was  remotely  controlled  using  the  LabView 
program  with  the  control  signals  being  sent  through  the  SCXI-1 124  module.  The  actual 
voltage  and  current  supplied  were  read  into  the  LabView  program  through  the  SCXI- 
1302  module  and  collected  by  the  slow  DAC.  The  high  voltage  was  placed  on  the  long 
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probe  making  it  an  ion  generator.  The  voltages  seleeted  for  eaeh  pressure  were  slightly 
lower  than  the  required  voltage  for  the  probe  to  are  to  the  wall  of  the  pressure  vessel.  Ion 
levels  were  then  measured  aeross  the  medium  probe  using  the  pieoammeter.  For  a  variety 
of  pressures,  the  ion  eurrent  was  examined.  Unfortunately,  the  extremely  low  eurrents 
generated  by  this  method  were  not  deteetable  by  the  pieoammeter.  Either  the  eurrent  was 
too  low  or  the  surrounding  noise  in  the  system  masked  any  measurable  eurrent  not  due  to 
the  eleetrie  field.  Sinee  this  investigation  was  not  eritieal  to  the  overall  understanding  of 
the  ion  sensor,  this  part  of  the  experiment  was  abandoned  in  favor  of  the  following 
method. 

III.IO  Transient  Ion  Current  Experiment 

Due  to  the  laek  of  results  from  the  steady  state  tests,  a  variant  of  that  method  was 
eondueted  by  investigating  the  unsteady  ease.  By  instantaneously  removing  the  high 
voltage  and  thereby  the  dominant  eleetrie  field,  the  remaining  eurrent  would  be  due  to  the 
ionization  of  the  gas  and  should  deeay  aeeording  to  the  pressure.  Again,  this  experiment 
made  use  of  a  high  voltage  power  supply  and  the  LabView  program.  By  elieking  on  a 
button  in  the  LabView,  the  Classman  power  supply  was  eommanded  to  0.0  V  and  the 
DACs  were  triggered  to  aequire  data.  Even  though,  theoretieally,  the  ion  deeay  rate  eould 
be  investigated  in  this  manner  the  Classman  power  supply  eould  not  instantaneously  drop 
the  voltage.  The  stored  energy  in  the  power  supply  drained  over  several  seeonds.  The 
resulting  ehange  in  eurrent  aeross  the  ion  probe  was  simply  due  to  the  ehange  in  the 
eleetrie  field.  Again,  the  physies  of  the  ion  deeay  eould  not  be  examined  due  to  the 
limitations  of  the  experimental  setup.  Although  a  different  eleetrieal  eireuit  eould  have 
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helped  alleviate  the  problem,  ion  generation  levels  were  still  expeeted  to  be  diffieult  to 
measure  without  major  experiment  redesign.  This  method  was  also  abandoned  in  favor  of 
the  eombustion  portion  of  the  experiments. 

III.ll  Combustion  Experiment 

The  main  effort  in  this  investigation  was  the  eombustion  experiment  where  the 
ion  eurrent  levels  generated  by  the  eombustion  of  methane  and  dry  air  could  be  directly 
compared  to  the  pressure. 

IILll.l  Initial  Conditions 

A  band  heater  temperature  of  50  degrees  Celsius  (°C)  created  the  buoyancy 
needed  for  mixing  of  the  methane  and  dry  air.  This  band  heater  temperature  corresponded 
to  an  internal  temperature  of  approximately  30  °C.  The  fuel-air  mixture  was  difficult  to 
ignite  at  1  atm  at  this  temperature.  A  higher  temperature  of  99  °C  increased  the  mixing 
within  the  vessel  and  eased  ignition  difficulties.  This  band  heater  setting  corresponded  to 
an  internal  temperature  of  45  to  65  °C.  Successive  tests  caused  the  heating  of  the  vessel 
and  increased  the  initial  temperature  despite  removal  and  injection  of  new  reactants. 

Some  variation  therefore  occurred  in  the  initial  temperature.  All  tests  began  with  the  band 
heater  setting  at  99  °C. 

IIL11.2  Test  Cases 

The  baseline  case  was  considered  the  test  with  an  initial  pressure  of  3  atm, 
equivalence  ratio  of  1 .0,  using  the  medium  probe  for  the  ion  sensor  with  a  +10  V  DC 
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bias,  a  band  heater  setting  of  99  °C,  and  the  pieoammeter  seale  set  to  20  pA.  The 
following  test  matriees  were  developed  based  on  searehing  out  from  this  baseline  test. 


Table  1.  Test  eases  for  medium  probe  at  +10  V  DC 

Equilvance  Ratio  ((p) 


Pressure  (atm) 

0.7 

1 

1.2 

1 

3/3 

3/1 

3/1 

2 

1/1 

1/0 

1/1 

3 

3/3 

3/3 

3/3 

5 

1/1 

3/3 

1/1 

7 

1/1 

Number  of  tests  at  High  Scale/Number  of  tests  at  low 

scale 


Table  2.  Test  eases  for  long  probe  at  +10  V  DC 

Equilvance  Ratio  (cp) 


Pressure  (atm) 

0.7 

1 

1.2 

1 

2 

3 

1/1 

2/1 

1/1 

5 

7 

Number  of  tests  at  High  Scale/Number  of  tests  at  low 
scale 
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Table  3.  Test  cases  for  short  probe  at  +10  V  DC 

Equilvance  Ratio  ((p) 


Pressure  (atm) 

0.7 

1 

1.2 

1 

2 

3 

1/1 

1/1 

1/1 

5 

1/0 

7 

Number  of  tests  at  High  Scaie/Number  of  tests  at  iow 
scaie 


Table  4.  Test  cases  at  various  voltages  at  3  atm  and  cp  =  1.0 

Probe 


Voitage  (DC) 

Short 

Medium 

Long 

5 

1/1 

10 

1/1 

3/3 

2/1 

20 

1/0 

40 

1/1 

1/0 

Number  of  tests  at  High  Scaie/Number  of  tests  at  iow 
scaie 
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111.11,3  Test  Procedure 


The  procedure  for  running  each  test  was  as  follows.  The  desired  test  case  was 
selected  and  the  fuel  and  air  procedure  was  used  to  fdl  the  vessel  with  the  appropriate 
partial  pressures  of  both  methane  and  dry  air.  The  power  supply  for  the  ignition  circuit 
relay  was  brought  to  12  V  DC  and  the  ignition  power  supply  was  raised  to  the  required 
12  V  DC.  After  ensuring  proper  conditions  with  the  continuous  mode  of  the  Lab  View 
program  the  mode  was  changed  to  data  acquisition.  After  confirming  the  program 
successfully  transferred  modes  and  waiting  a  couple  seconds,  the  ignition  switch  was 
pressed  causing  a  spark  that  ignited  the  mixture.  Because  a  small  but  noticeable  amount 
of  noise  was  generated  simply  by  depressing  the  switch,  the  momentary  switch  was  not 
released  until  after  1  second.  This  ensured  any  noise  generated  by  the  switch  would  not 
influence  the  results.  The  timing  of  the  ignition  also  coincided  with  no  current  following 
through  the  band  heater.  This  was  accomplished  by  visually  watching  the  indictor  on  the 
band  heater  controller.  The  data  files  were  then  saved  to  the  hard  drive.  The  combustion 
products  were  removed  following  the  previously  described  procedure  for  the  fuel  and  dry 
air. 
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IV.  Raw  Data 


IV.  1  Chapter  Overview 

This  chapter  describes  the  error  in  measuring  the  physieal  values  and  shows  the 
raw  data  before  any  data  reduetion.  Interesting  features  of  the  raw  data  are  highlighted. 

IV.2  Measurement  Error 

For  each  measurement  the  error  was  ealeulated  from  the  aeeuraeies  of  eaeh  deviee 
involved  in  the  measurement.  The  time  error  in  the  eurrent  measurements  was  one  half  of 
the  time  between  eaeh  sample  beeause  of  the  aeeurate  timing  of  the  hardware  trigger  of 
the  fast  DAC  as  previously  diseussed.  The  time  error  for  the  slow  DAC,  also  as 
previously  diseussed,  was  not  only  half  the  sample  time  but  ineluded  possible  errors  from 
the  eonditional  triggers.  Assuming  the  errors  imposed  by  eonditional  trigger  are 
negligible,  the  time  errors  were  0.1  ms  and  10  ps  for  the  pressure  and  eurrent 
respeetively. 

The  root  of  the  sum  of  the  squares  (RSS)  method,  as  deseribed  by  Wheeler  and 
Ganji  [24],  provided  a  suitable  error  estimate  for  the  pressure  and  eurrent  measurements. 
The  error  in  the  pressure  measurements  was  0.0513  atm  for  initial  pressures  of  latm  to  3 
atm  and  an  error  of  0.171  atm  for  initial  pressures  of  5  atm  and  7  atm.  The  error  in  the 
eurrent  measurement  was  2.79%  for  the  20  pA  seale. 
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For  each  figure,  error  bars  were  added  at  several  points  to  indicate  the  error  levels, 
how  the  variations  in  each  test  compared  to  the  error,  and  how  the  signal  noise  compared 
to  the  error  levels. 

IV,3  Results  of  Raw  Data 

The  baseline  case  with  an  equivalence  ratio  of  1 .0  and  initial  pressure  of  3  atm 
was  repeated  3  times  as  shown  in  Figure  1 1 .  Unless  otherwise  stated,  each  case  used  the 
medium  probe  energized  with  10  V  DC  with  a  band  heater  setting  of  99°  C  and  a  scale  of 
20  pA  on  the  picoammeter.  Time  of  0.0  seconds  corresponds  to  depressing  the  ignition 
switch  and  not  the  time  the  spark  actually  occurred. 
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-0.1  0  0.1  0.2  0.3  0.4  0.5 

Time  (s) 

Figure  1 1 .  Runs  at  the  baseline  ease  of  cp  =  1 .0  at  an  initial  pressure  of  3  atm 


The  large  noise  just  after  0.0  seeonds  is  due  to  interferenee  and  effeets  of  the  ignition 
spark.  The  square  and  eirele  objeets  are  indietors  of  where  errors  were  eomputed  and  not 
the  aetual  size  of  the  error.  Error  bars  within  the  objeets  show  the  aetual  size  of  the  error. 
The  measurement  errors  in  Figure  1 1  are  diffieult  to  diseern  due  to  their  relatively  small 
size  eompared  to  aetual  measurement  levels.  Figure  1 1  overall  shows  good  repeatability 
of  the  measurements. 
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Pressure  (atm) 


The  impact  of  changing  the  ion  probe  voltage  from  10  V  DC  to  40  V  DC  and  5  V 


DC  is  shown  in  Figure  12. 
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Figure  12.  Comparison  of  5V,  lOV,  &  40V  across  the  ion  probe  at 
an  initial  pressure  of  3  atm  and  (p  =  1 .0 


The  5V  case  shows  reduced  current  levels  but  similar  trends  to  the  lOV  case.  The  polarity 
of  the  current  switches  in  the  40V  case  and  the  signal  appears  to  become  saturated. 

Again,  error  bars  are  shown  at  select  points  for  comparison  purposes. 


56 


Pressure  (atm) 


At  an  initial  pressure  of  1  atm,  the  effeets  of  ehanging  the  equivalenee  ratio  are 
shown  in  Figure  13. 
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Figure  13.  Effeets  of  (p  =  0.7,  1.0,  &  1.2  at  initial  pressure  of  1  atm 


Both  the  (p  =  1.0  and  1.2  eases  show  similar  measurement  levels  and  trends.  The  (p  =  0.7 
case  shows  a  time  delay  in  the  pressure  peak  and  lower  peak  level.  The  current  in  this 
case  does  not  display  a  sharp  spike  but  a  low  peak  of  approximately  2  pA. 
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Pressure  (atm) 


At  an  initial  pressure  of  3  atm,  the  effeets  of  ehanging  the  equivalenee  ratio  are 


shown  in  Figure  14. 
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Figure  14.  Effeets  of  (p  =  0.7,  1.0,  &  1.2  at  initial  pressure  of  3  atm 

Again,  both  the  (p  =  1.0  and  1.2  eases  show  similar  measurement  levels  and  trends.  The  (p 
=  0.7  ease  shows  a  time  delay  in  the  pressure  peak  and  lower  peak  level.  The  eurrent  in 
this  ease  does  not  display  a  sharp  spike  but  a  low  peak  of  approximately  1  pA. 
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Pressure  (atm) 


At  an  initial  pressure  of  3  atm,  the  effeets  of  ehanging  the  equivalenee  ratio  are 


shown  in  Figure  15. 
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Figure  15.  Effeets  of  cp  =  0.7,  1.0,  &  1.2  at  initial  pressure  of  5  atm 


As  previously  shown  in  Figures  3  and  4,  the  (p  =  1.0  and  1.2  eases  show  similar 
measurement  levels  and  trends.  The  cp  =  0.7  ease  shows  a  time  delay  in  the  pressure  peak 
and  lower  peak  level.  The  eurrent  is  this  ease  does  not  display  a  sharp  spike  but  a  low 
peak  of  approximately  0.5  pA. 
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Pressure  (atm) 


The  changes  in  pressure  and  ion  current  for  initial  pressures  of  1,  3,  &  5  aims  are 
shown  in  Figure  16. 
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-0.1  0  0.1  0.2  0.3  0.4  0.5 

Time  (s) 

Figure  16.  Effect  of  initial  pressures  of  1,  3,  &  5  atm  at  cp  =  1.0 

As  seen  in  Figure  6,  increasing  initial  pressure  results  in  lower  ion  current  peaks.  The 
trends  of  the  ion  current  and  instantaneous  pressure  for  each  initial  pressure  are  similar 
however  the  timing  of  the  pressure  peak  shows  a  slight  time  delay  at  higher  initial 
pressures. 
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Pressure  (atm) 


The  area  of  interest  in  Figure  6  is  expanded  and  more  elearly  shown  in  Figure  17. 


Figure  17.  Zoomed  in  view  of  Figure  16,  effeets  of  initial  pressure  of  1,  3,  &  5 

atm  at  (p  =  1 .0 

In  Figure  17,  the  ion  eurrent  also  shows  a  slight  time  delay  for  higher  pressures. 

The  effeets  of  changing  the  probe  length  for  the  baseline  case  of  3  atm  at  cp  =  1.0 
are  shown  in  Figure  18. 
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Figure  18.  Effects  of  probe  length  at  baseline  case  of  cp  =  1.0  at  3  atm 

The  large  narrow  spike  near  0.0  seconds  was  due  to  the  interference  of  the  ignition  spark. 
The  baseline  medium  probe  produced  the  highest  ion  current  peak  followed  by  the  long 
probe  with  the  short  probe  producing  the  smallest  ion  current  peak.  The  short  probe 
current  peak  had  a  slight  time  delay  compared  to  the  medium  probe  due  to  the  further 
distance  from  the  ignition  spark.  Likewise,  the  long  probe  showed  a  greater  time  delay 
since  the  1 0  cm  extension  placed  it  even  further  away  from  the  source  of  the  ignition. 
Figure  19  shows  an  exploded  view  of  the  features  of  interest  in  Figure  18. 
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Figure  19.  Zoomed  in  view  of  Figure  18  more  clearly  showing  the  effects  of 
probe  length  at  the  baseline  case  of  cp  =  1.0  at  3  atm 


The  pressure  measurements  in  Figure  9  show  relatively  small  variations  while  the  ion 
current  peak  shows  large  differences  in  both  peak  level  and  peak  timing.  In  Figure  19, 
the  time  delay  between  the  ion  current  peaks  is  more  clearly  evident. 
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Pressure  (atm) 


V.  Data  Reduction  and  Discussion 


V.l  Chapter  Overview 

In  this  chapter  the  data  is  analyzed  and  eompared  against  predietions  of  Chapter 
II.  Features  of  the  raw  data  are  examined  and  eompared  against  expeetations.  The  data  is 
then  redueed  by  several  methods  ineluding  averaging  the  results  over  similar  eases, 
eomputing  the  derivative  of  the  ion  eurrent  signal,  and  filtering  out  noise.  Features  of  the 
redueed  data  are  also  examined  and  eompared  against  expeetations. 

V,2  Analysis  of  Raw  Data 

The  ion  eurrent  from  the  raw  data  shown  in  Chapter  IV  displays  initial  noise  due 
to  the  ignition  spark,  a  sharp  peak  a  short  time  later  due  to  the  chemi-ionization  as  the 
flame  front  passes  the  probe,  and  a  smaller  seeondary  peak  after  the  sharp  peak.  The 
strength  of  the  seeondary  peak  is  highest  for  the  3  atm  ease  as  shown  in  Figure  16.  When 
eompared  to  the  sharp  ehemi-ionization  spike,  the  shape  and  strength  of  the  seeondary 
peak  is  irregular  even  for  repeated  eases  as  shown  in  Figure  11.  This  seeondary  peak 
could  be  the  result  of  NO  production  similar  to  the  IC  engines.  Sinee  both  temperature 
and  pressure  are  mueh  lower  than  in  IC  engines,  the  NO  would  not  be  the  dominant 
souree  of  ionization  in  the  seeondary  peak.  Non-uniformities  and  a  three  dimensional 
flame  front  are  more  likely  the  eause  of  the  secondary  peaks.  Due  to  the  surfaee  area  of 
the  medium  probe,  the  flame  front  may  pass  near  the  sensor  multiple  times  producing 
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additional  ionization.  The  variability  and  structure  may  provide  useful  information  for 
PDE  but  that  investigation  is  beyond  the  scope  of  this  thesis. 

Equivalence  ratios  of  1.0  and  1.2  show  similar  signal  levels  and  trends.  The  0.7 
equivalence  ratio,  however,  shows  little  resemblance  to  the  other  two  ratios.  The  pressure 
peak  is  significantly  lower  and  the  ion  current  spike,  the  dominant  feature  of  the  other 
two  ratios,  is  nonexistent.  This  trend  of  lower  ion  currents  with  low  equivalence  ratios  is 
expected.  Equation  (1)  implies  the  ionization  production  strongly  depends  upon  the 
concentration  of  the  fuel.  Eurther,  the  recombination  of  ions,  as  shown  in  Eq  (5),  strongly 
depends  upon  the  oxygen  concentration.  The  complete  absence  of  the  spike,  although  not 
expected,  agrees  with  the  trend  of  lower  ion  currents  at  lower  equivalence  ratios.  In 
addition,  the  only  peak  in  the  0.7  equivalence  ratio  cases  correspond,  although  a  lower 
level,  to  the  secondary  peaks  shown  in  cases  of  1.0  and  1.2  equivalence  ratio.  This  also 
agrees  with  the  expected  trend. 

Since  the  Omega  pressure  transducer  had  a  stand  off  distance  of  approximately  10 
cm  from  the  lid  of  the  pressure  vessel,  this  could  impact  the  pressure  measurements.  Any 
stand  off  distance  would  influence  the  measurements  by  smoothing  any  pressure  spikes  in 
the  vessel.  Because  of  the  deflagration  in  the  constant  volume  process,  pressure  spikes 
are  not  expected.  The  PDE,  however,  does  have  the  von  Neumann  pressure  spike  as 
previously  described  by  the  ZND  wave  structure.  Eor  the  purposes  of  this  experiment  any 
smoothing  of  the  pressure  measurement  is  considered  negligible.  The  stand  off  distance 
still  can  introduce  a  time  delay  in  the  pressure  signal.  At  an  approximate  initial  gas 
temperature  of  50  °C,  the  mixture  has  a  speed  of  sound  of  roughly  360  m/s.  A  stand-off 
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distance  of  10  cm  therefore  results  in  a  time  delay  about  0.028  ms.  This  small  delay  also 
can  be  considered  negligible. 


V.3  Error  Analysis 

In  addition  to  the  measurement  errors  calculated  in  Chapter  IV,  the  error  of  the 
equivalence  ratio  was  also  computed  using  the  RSS  method.  The  resulting  error  for  each 
initial  pressure  and  equivalence  ratio  is  shown  in  Table  5. 


Table  5.  Percent  Error  in  Equivalence  Ratio 


Equivalence  Ratio 

Pressure  (atm) 

0.7 

1.0 

1.2 

1 

11.2% 

11.9% 

12.3  % 

2 

5.9  % 

3 

3.7  % 

4.0  % 

4.1  % 

5 

2.2  % 

2.4  % 

2.5  % 

7 

1.7% 

The  errors  at  1  atm  are  greater  due  to  the  limited  accuracy  of  the  15  psia  Endecvo 
transducer  for  both  the  vacuum  pressure  and  the  partial  pressure  of  the  fuel.  At  higher 
total  pressures,  the  relative  error  in  the  fuel  partial  pressure  decreases. 
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Averaged  Current  (MicroAmps) 


In  analyzing  the  ion  current  the  derivative  of  the  signal  was  calculated  using  a 
fourth  order  accurate  scheme.  The  resulting  error  in  the  derivative  was  3.94  %  compared 


to  2.79  %  for  the  ion  current  signal. 

V,4  Reduced  Data  as  a  Function  of  Time 

Because  of  the  variability  of  the  ion  current,  the  average  of  three  runs  at  one 
condition  was  calculated  as  shown  in  Figure  20. 
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Figure  20.  Averaged  Current  and  Pressure  for  initial  pressures  of  1,  3,  &  5  atm 

at  (p  =1.0 
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Averaged  Pressure  tatmt 


As  in  Chapter  IV,  error  bars  are  ineluded  in  Figure  20  at  seleet  points.  The  time  seale  was 
narrowed  to  allow  better  resolution  into  the  eurrent  spike.  Compared  to  the  signal  levels, 
the  error  in  both  the  eurrent  and  pressure  measurements  are  relatively  small.  The  initial 
noise  in  the  ion  eurrent  is  due  to  the  ignition  spark.  When  averaged  over  the  three  runs, 
the  time  differenees  between  the  signals  ereate  a  smoothing  of  the  eurrent  spike.  Despite 
this  effeet,  the  same  trends  are  seen  with  the  1  atm  ease  showing  the  strongest  eurrent 
spike  and  the  5  atm  ease  showing  the  weakest  eurrent  spike. 

The  pressure  and  eurrent  ratio  model  developed  by  Saitzkoff  et  al.  was  applied  to 
the  baseline  test  condition  to  evaluate  the  suitability  of  the  model  to  these  test  conditions. 
After  normalizing  the  measured  pressure,  the  predicted  current  was  plotted  in  Figure  21 
and  compared  to  the  measured  current  after  normalization. 

The  model  was  developed  for  use  in  IC  engines  where  the  maximum  pressure  and 
combustion  occur  when  the  piston  is  near  the  top  of  cylinder.  The  constant  volume 
combustion  process  is  inherently  different.  In  the  constant  volume  process  the  flame  front 
passes  the  sensor  well  before  the  maximum  pressure  in  the  vessel  is  produced.  This  effect 
is  clearly  shown  in  Figure  21  by  showing  the  measured  current  spike  well  before  the 
predicted  values.  Based  on  these  results,  the  Saitzkoff  et  al.  model  is  not  suitable  for 
pressure  predictions  in  the  constant  volume  process. 
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Normalized  Value 


Figure  21.  Saitzkoff  et  al.  model  applied  to  this  experiment  at  baseline 
condition  of  3  atm  initial  pressure  and  (p  =  1 .0 


The  rate  of  change  or  derivative  of  the  ion  current  was  predicted  in  chapter  11  to 
be  a  function  of  the  pressure  squared.  The  derivative  was  calculated  using  a  4*''  order 
accurate  scheme: 


-M/+2 

Uh 
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where  h  =  At.  and  i  is  the  point  the  about  whieh  the  derivative  is  eomputed.  By  taking  the 
derivative,  the  noise  was  amplified  eompared  to  the  original  ion  eurrent  as  shown  in 
Figure  22. 


Figure  22.  Current  and  derivative  at  baseline  ease 

The  resulting  noisy  derivative  makes  diseeming  useful  information  diffieult  at  best.  A 
limited  time  seale  permits  better  resolution  into  the  region  of  interest. 

In  order  to  extraet  information,  a  low-pass  finite  impulse  response  filter  was 
applied  using  the  remez  funetion  in  MatLab.  This  filter  was  designed  by  using  a 
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transition  zone  from  500  Hz  to  1  kHz.  In  this  manner  unwanted  noise  above  1  kHz  is 
eliminated  but  important  oharaeteristies  of  the  signal  remain  as  shown  in  Figure  23. 
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Figure  23.  Filtered  derivative  and  eurrent  at  baseline  ease 

Again  the  time  seale  is  redueed  to  resolve  the  rate  of  ehange  from  the  maximum  eurrent. 
By  filtering  the  derivative,  the  noise  is  redueed  but  still  shows  some  oseillatory  behavior. 
The  error  bars  show  that  although  the  error  inereases  by  taking  the  derivative,  the  error  is 
still  relatively  small  eompared  to  the  overall  derivative  level. 
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The  repeatability  of  the  filtered  derivative  is  shown  in  Figure  24.  Again  some  time 
variation  oeeurs  between  eaeh  of  the  three  runs.  Peak  levels  also  very  between  runs  but 
similar  behavior  and  patterns  indict  that  signal  is  not  random  but  contains  useful 
reproducible  information. 


K) 
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Figure  24.  Repeatability  of  filtered  derivative  for  3  runs  at  baseline  condition 
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Derivative  of  Current  (MicroAmp/s) 


How  the  derivative  is  affected  by  pressure  is  investigated  in  Figure  25.  The  area  of 
interest  in  expanded  in  Figure  26. 
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Figure  25.  Filtered  derivative  eompared  for  several  initial  pressures  at  (p  =  1.0 
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Figure  26.  Expanded  view  of  Figure  25  showing  pressure  affeets  on  the 

derivative 


The  maximum  positive  level  of  the  derivative  clearly  decreases  at  higher  pressures.  The 
minimum  value  of  the  derivative  is  not  as  clear  but  tends  to  increase  with  pressure.  That 
is,  the  rate  of  ion  current  decay  from  the  peak  decreases  with  pressure. 

V.5  Influence  of  Equivalence  Ratio 

Because  the  ion  current  depends  upon  the  chemi-ionization  process  during 
combustion,  the  equivalence  ratio  is  a  critical  factor  for  the  ionization  levels.  As 
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previously  discussed,  the  ion  current  was  expected  to  decrease  with  lower  equivalence 
ratios.  The  raw  data  presented  in  Chapter  IV  agreed  with  the  trend  to  the  extent  that  an 
ionization  spike  was  nonexistent.  By  comparing  the  peak  ion  currents,  this  trend  can  be 
better  understood  as  shown  in  Figure  27. 


Figure  27.  Effect  of  equivalence  ratio  on  maximum  ion  current  for  initial 

pressures  of  1,  3,  &  5  atm 

Because  of  the  high  relative  error  in  the  fuel  partial  pressure  at  1  atm,  the  error  in 
equivalence  ratio  is  substantial.  Despite  this  uncertainty,  the  ion  current  decreases  sharply 
with  lower  equivalence  ratios.  At  higher  equivalence  ratios,  the  ion  current  decreases 
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slightly  for  the  3  atm  and  5  atm  cases  but  increases  for  the  1  atm  case.  Since  ion 
production  requires  both  fuel  and  oxygen  from  equation  1 ,  an  ionization  peak  is  expected 
near  or  above  stoichiometric  conditions.  The  limited  number  of  tests  run  at  various 
equivalence  ratios  provides  little  resolution  into  the  exact  location  of  this  peak  and  how 
the  location  is  affected  by  pressure. 

The  derivative  of  the  current  as  a  function  of  equivalence  ratio  can  be  examined  in 
a  similar  fashion  as  shown  in  Figure  28. 


Figure  28.  Peak  decay  rate  of  current  as  a  function  of  equivalence  ratio  for 
initial  pressures  of  1,  3,  &  5  atm 
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The  peak  deeay  rate  or  minimum  negative  derivative  of  the  eurrent  oeeurs  at 
stoiehiometrie  eonditions  for  the  3  atm  and  5  atm  eases  but  at  an  equivalenee  ratio  of  1 .2 
for  the  1  atm  ease.  Compared  to  the  maximum  eurrent,  the  peak  deeay  rate  shows  the 
similar  trends.  That  is,  where  the  ion  eurrent  was  a  maximum  at  stoiehiometrie 
eonditions,  the  peak  deeay  rate  is  also  greatest  at  stoiehiometrie  eonditions. 

The  noise  inherent  in  the  derivative  prevents  useful  examination  of  individual 
points.  An  average  deeay  rate  over  a  seleeted  range  also  has  limitations  and  is  somewhat 
arbitrary.  The  peak  ion  eurrent  provides  a  elear  starting  point  for  the  range,  but  a  finding 
a  suitable  end  point  is  challenging.  Selecting  the  end  point  to  be  when  the  derivative 
equals  zero  creates  easy  calculations  but  the  large  time  variance  between  test  cases 
directly  affects  the  result.  End  points  based  on  other  criteria  are  also  highly  variable 
because  of  both  the  noise  and  time  variation  in  the  signal.  The  parameter  chosen  to  best 
quantify  the  derivative  was  the  peak  decay  rate  because  of  the  relative  consistency  over 
several  test  runs. 

V,6  Results  as  a  Function  of  Pressure 

Ultimately,  the  objective  is  to  calculate  the  pressure  by  measuring  the  ion  current. 
The  theory  of  Chapter  II  states  that  the  rate  of  change  of  the  ion  current,  or  decay  rate, 
should  be  a  function  of  pressure.  The  raw  results  of  Chapter  IV  show  that  pressure  has  a 
strong  influence  on  the  peak  ion  current.  The  peak  ion  current  averaged  over  three  runs  is 
compared  to  the  pressure  and  shown  in  Figure  29. 
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Figure  29.  Maximum  current  averaged  over  three  test  runs  plotted  against  the 
instantaneous  pressure  at  cp  =  1 .0 


The  linear  trend  line  in  Figure  29,  computed  by  the  least  squares  method,  shows  a  strong 
correlation  with  an  R-squared  value  of  0.9982. 

The  computed  linear  trend  is  then  compared  to  the  individual  cases  as  shown  in 
Figure  30. 
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Figure  30.  Maximum  current  for  individual  runs  compared  to  pressure.  The 
trend  line  shows  the  correlation  computed  from  averaged  cases 

When  compared  against  the  individual  test  runs,  the  linear  trend  line  is  far  less  accurate. 

The  same  approach  is  repeated  for  the  derivative  of  the  ion  current.  The  peak 
decay  rate  averaged  over  three  test  runs  is  compared  to  the  pressure  as  shown  in  Figure 
31. 
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Figure  31.  Peak  deeay  rate  averaged  over  three  test  runs  and  eompared  against 
instantaneous  pressure  for  oases  of  (p  =  1.0 

Again,  a  linear  trend  line  oomputed  with  the  least  squares  method  shows  a  good 
oorrelation  between  the  pressure  and  the  averaged  peak  deeay  rate  with  an  R-squared 
value  of  0.9997.  Using  this  relationship  to  determine  pressure,  the  resulting  linear  model 
is: 


R  =  0.0043964 


^  mean  ) 

dt 


+  9.491 


(22) 
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is  the  averaged  peak  deeay  rate  in  pA/s  and  P  is  the  predieted  pressure  in 


dl 
dt 

atm.  This  linear  model  is  only  applieable  to  a  pressure  range  from  1  atm  to  5  atm  at 
stoichiometrie  conditions.  The  usefulness  of  the  model  is  also  limited  because  the  mean 
peak  decay  rate  is  required.  This  value  is  rarely  available  since  actual  conditions  can  vary 
for  each  measurement. 

A  more  useful  model  can  be  devised  by  examining  individual  cases  and  not  the 
mean  peak  decay  rate.  Again  a  linear  correlation  trend  line  is  computed  and  the  result,  as 
shown  in  Figure  32,  is  similar  to  the  previous  correlation  but  the  individual  cases  show 
greater  variability. 
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Figure  32.  Peak  deeay  rate  for  individual  eases  eompared  to  the  instantaneous 
pressure.  A  new  eorrelation  is  eomputed. 

The  linear  eurve  fit  with  an  R-squared  value  of  0.7455  is  not  as  aeeurate  as  the  mean 
ease.  The  variability  of  the  peak  deeay  rate  with  eaeh  test  reduees  the  aeeuraey  of  the 
linear  eurve  lit.  A  more  useful  linear  model  based  on  these  individual  eases  is: 

R  =  0.0047515— +  10.084  (23) 

dt 


82 


where  —  is  the  peak  decay  rate  in  pA/s  and  P  is  the  predicted  pressure  in  atm.  This 

linear  model  is  only  applicable  to  a  pressure  range  from  1  atm  to  7  atm  at  stoichiometric 
conditions. 

The  root  mean  square  error  (RMSE)  is  one  useful  metric  to  describe  the  error  of 
the  linear  model.  The  definition  of  RMSE  is: 


RMSF=  -X(x,-ay  (24) 

V «  i=i 

where  n  is  the  number  of  samples,  Xi  is  the  value  at  sample  i  and  a  is  the  mean.  Elsing  the 
statistics  software  package  IMP,  the  RMSE  calculated  for  the  linear  model  is  1.0635  atm. 
Although  the  accuracy  of  this  linear  model  is  less  than  desired,  the  model  does  show  a 
strong,  predictable  relationship  between  the  peak  decay  rate  and  pressure. 

The  linear  model  includes  an  intercept  term  that  was  not  predicted  by  equations 
19  &  20.  Since  the  intercept  term  is  on  the  same  order  as  the  slope  term  it  can  not  be 
neglected.  This  term  may  be  due  to  increasing  ion  recombination  rates  at  higher  pressures 
preventing  the  chemi-ionization  from  producing  a  super-equilibrium  ion  concentration. 
Near  10  atm,  the  ion  recombination  rate  may  match  or  exceed  the  ion  production  rate. 
This  agrees  with  the  expectation  that  the  ion  decay  rate  is  a  function  of  the  square  of 
pressure.  This  may  limit  the  range  of  pressures  that  can  be  measured  by  this  method. 
Nonetheless,  this  range  is  sufficient  for  capturing  important  pressures  in  the  PDE. 
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V.7  Probe  Resolution 


The  ion  probe  speeial  resolution  is  on  the  same  order  of  the  diameter  of  the  probe, 
approximately  0.3  em.  The  time  resolution  of  the  probe  can  be  approximated  by; 


t  = 


(25) 


where  t  is  the  time,  Vavg  is  the  speed  of  sound  (360  m/s),  m  is  the  mass  of  a  molecule 
(4.84  X  10'^^  kg  for  CHO^),  q  is  the  elementary  charge  (1.60  X  10'^^  C),  and  E  is  the 
applied  electric  field  of  1 .0  kV/m.  The  resulting  response  time  is  on  the  order  of  100 
nanoseconds  or  10  MegaHertz  (MHz). 

The  time  and  spatial  resolution  of  the  probe  is  adequate  to  observe  the  detonation 
wave  structure  and  the  following  expansion  zone  in  the  PDE  but  not  the  von  Neumann 
pressure  spike. 
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VI.  Conclusions 


VLl  Conclusions 

The  objective  stated  in  Chapter  I  is  to  develop  a  new  pressure  sensor  for  the  PDE 
based  on  work  already  accomplished  in  IC  engines  and  the  PDE  itself  Chapter  II 
predicts  that  the  rate  of  change  of  the  ion  current  decay  rate  should  be  a  linear  function  of 
pressure.  Inversely,  pressure  should  be  a  linear  function  of  the  ion  current  decay  rate.  The 
results  of  Chapters  IV  and  V  show  that  both  the  ion  current  and  the  peak  decay  rate  are  a 
linear  function  of  pressure.  In  other  words,  pressure  is  a  linear  function  of  the  ion  current 
and  the  peak  decay  rate.  The  theoretical  prediction  is  generally  correct  in  predicting  a 
linear  relationship  despite  variations  in  the  ion  current  decay  rate.  The  developed  linear 
model  in  equation  23  provides  an  initial  useful  model  for  determining  pressure  although 
the  accuracy  needs  to  be  improved.  The  intercept  term  in  equation  23  may  limit  the  useful 
range  of  this  method. 

The  variations  can  be  explained  by  several  factors  including  local  changes  in 
equivalence  ratio  and  non-uniform  flame  fronts.  The  derivation  also  did  not  account  for 
additional  forces  and  phenomena  that  can  affect  both  the  ion  decay  rate  and  the  drift 
velocity.  The  instrumentation  in  the  experimental  setup,  particularly  the  picoammeter, 
may  not  have  been  fast  enough  to  correctly  capture  the  transient  events. 

Interestingly,  both  the  peak  decay  rate  and  the  actual  ion  current  showed  a  linear 
relationship  with  pressure.  The  IC  engines  only  used  the  ion  current  to  determine  the 
pressure.  As  previously  discussed,  large  variations  in  current  can  occur  due  to  a  number 
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of  factors  especially  impurities.  The  derivative  or  peak  decay  rate  method  inherently 
inereases  the  noise  in  the  signal  but  offers  a  means  of  addressing  problems  with 
variations  in  the  aetual  ion  eurrent.  Although  improved  aeeuraey  is  required,  this 
investigation  eonceptually  shows  the  spark  plug  can  provide  an  alternate  method  to 
measuring  pressure  in  the  PDE. 

VL2  Future  Work 

Future  work  will  need  a  more  detailed  examination  of  the  phenomena  ereating  the 
relationship  between  pressure  and  the  peak  deeay  rate.  Additional  testing  will  be  required 
to  address  other  eonsiderations  sueh  as  different  hydroearbon  fuels  and  impurities  in  the 
fuel.  A  larger  set  of  test  cases  would  help  improve  the  usefulness  of  this  method  and 
aeouraey  of  the  model.  Ultimately,  the  sensor  will  need  to  be  compared  to  eurrent 
pressure  sensing  methods  within  the  PDE. 

VI,3  Future  Experimental  Setup  Recommendations 

While  the  test  setup  was  adequate  for  this  initial  investigation,  subsequent  work 
should  improve  this  testing  approaeh.  Inereasing  the  sample  rate  of  the  picoammeter 
would  improve  eonfidenee  that  the  transient  events  are  eorreetly  eaptured.  This 
improvement  will  also  be  needed  for  the  faster  processing  in  the  detonation  cyele 
eompared  to  the  eonstant  volume  proeess  investigated  in  this  experiment.  Also,  to  ease  in 
data  eapture,  a  DAC  capable  of  hardware  triggering  should  replace  the  E-series  card.  This 
will  release  eomputer  resourees  for  better  time  resolution  and  improve  time  aeeuraey  not 
only  on  that  eard  but  with  the  faster  NI  61 10  card  as  well.  A  faster  computer  is  always 
desired  but  in  this  testing  the  limited  256  MB  of  memory  was  the  biggest  detriment. 
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The  K-type  thermocouple  response  time  was  too  slow  to  provide  much  useful 
thermodynamic  information  other  than  initial  temperature  in  the  combustion  bomb.  A 
thermocouple  with  a  smaller  bead  size  would  be  much  more  useful. 

The  large  errors  in  equivalence  ratio  at  low  pressure  could  easily  be  reduced  by 
using  additional  pressure  transducers  to  more  accurately  determine  the  partial  and  total 
pressures.  A  5  psia  transducer  would  be  effective  for  measuring  the  vacuum  and  fuel 
partial  pressures  at  low  initial  pressures.  A  50  psia  would  be  useful  in  reducing  the  error 
in  the  total  pressure  for  the  2  atm  and  3  atm  test  cases. 

The  variation  in  probe  response  with  voltage  was  also  not  well  understood. 
Saturation  appeared  to  occur  but  the  exact  cause  was  not  determined.  Further 
investigation  into  the  cause  would  also  be  useful  in  development  of  this  method. 

The  poor  wiring  in  the  lab  made  creating  a  true  ground  impossible.  Although 
efforts  were  taken  to  reduce  noise  as  much  as  feasible,  the  experiment  was  still 
susceptible  to  noise  from  many  sources. 

Testing  for  a  larger  pressure  range  could  also  aid  understanding  and  performance 
of  the  ion  sensor.  A  larger  number  of  cases  at  one  condition  would  improve  statistical 
confidence  in  the  measurement. 

Although  some  effort  was  made  in  selecting  an  appropriate  filter  for  the  ion 
current  derivative,  this  method  was  not  optimized  nor  was  the  effect  of  the  filter  fully 
taken  into  account  during  the  data  reduction. 

The  developed  linear  model  works  fairly  well.  A  better  statistical  analysis  with 
more  data  points  could  reveal  a  more  accurate  model  as  well  as  provided  more  insight 
into  the  dominant  physical  phenomena  especially  the  intercept  term. 
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Appendix  A:  LabView  Program 


LabView  version  7.0  provided  a  quick  and  easy  software  package  to  collect  the 
data.  LabView  has  a  large  set  of  drivers  to  work  with  a  variety  of  instruments  and  data 
acquisition  cards.  Initially,  writing  a  LabView  code  was  expected  to  be  fairly  straight 
forward  because  of  all  of  its  capabilities.  Unfortunately,  no  pre-built  Virtual  Instruments 
(VI)s  existed  that  performed  all  the  required  functions  so  a  new  VI  had  to  be  coded.  The 
user  friendly  Vis  built  into  LabView  were  too  slow  and  limited  to  be  useful  for  this 
experiment.  A  combination  of  higher  level  Vis  was  required  to  create  the  object  oriented 
code.  Newer  drivers  (LabView  mx)  were  not  yet  written  for  the  fast  data  acquisition  card 
so  the  older  traditional  driver  and  Vis  had  to  be  used.  The  lack  of  hardware  triggering  in 
the  E-series  card  required  an  elaborate  VI  intensive  conditional  triggering  scheme  as 
mentioned  earlier.  These  complications  turned  the  code  into  a  substantial  effort. 

The  coded  program  operated  in  two  different  modes.  The  first  mode  allowed  for 
ensuring  proper  conditions  before  beginning  the  test.  This  mainly  used  continuous  data 
acquisition  Vis  running  at  low  data  capture  rates  and  displayed  the  results  on  the 
computer  screen  in  real  time  as  shown  in  Figure  33.  To  help  determine  sources  of  noise 
during  the  experimental  setup,  a  EFT  VI  was  used  to  show  the  current  from  the 
picoammeter  in  the  frequency  domain  as  shown  in  Figure  34.  A  simple  digital  band-gap 
filter  VI  was  also  tested  but  not  used  for  data  capture  to  reduce  noise  from  the  dominant 
60  Hz  sources.  The  detailed  block  diagram  for  the  first  mode  is  shown  in  Figure  35  and 
Figure  36. 

The  second  program  mode  used  a  hardware  trigger  VI  for  the  fast  DAC  and  a 
conditional  trigger  VI  for  the  E-series  DAC.  These  Vis  captured  a  set  number  of  data 
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points  before  and  after  the  trigger.  After  aequiring  all  the  data,  results  were  saved  to  a  file 
and  displayed  on  the  screen  as  shown  in  Figure  37.  Any  errors  the  program  encountered 
were  also  displayed.  A  time  limit  was  also  used  for  both  these  Vis  to  make  sure  the 
DACs  did  not  go  into  a  continuous  loop  looking  for  the  trigger  condition.  The  detailed 
block  diagram  for  the  second  mode  is  shown  in  Figures  38  and  Figure  39. 

The  program  was  switched  from  the  first  mode  to  the  second  mode  by  clicking  on 
a  Boolean  button  that  moved  the  operation  from  one  case  structure  to  the  other. 
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Figure  33.  First  program  mode  front  panel 
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Figure  34.  First  mode  front  panel  showing  ion  current  in  frequency  domain 


Figure  35.  Top  section  of  block  diagram  for  first  mode 
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Figure  36.  Bottom  section  of  block  diagram  for  first  mode 


file  Edit  Operate  lools  growse  Window  Help 


[iT]  I  13pt  Application  Font  •.  1 1  ||  ||  j  |  j 


auu.uu  1- 

Taking  Data? 


J 


J 


lal  Temp  HU  Vdtage 
156  0 

ob*  (V)  HU  Current 


Enable  Ualue 
.0195e" 


buffer  £t2e 

SilOOOOO  I 

number  of  scar 
to  acquire 

SlllQQO  I 


gllOOOO.OO  I 

pretrigger 


nurrtter  cf  sc 
to  acquire 
glllOOOO 


gaoOOOO.QOl 


.IJ 

1 

Sisiatll  ®  SS  ©  1 1  '■•“EW 

ll^EKpTb.vi  Front  Panel  I^Documentl  -  Mcrosoft  UU...  | 

|S.<f  *.a%  16.26 

Figure  37.  Front  panel  for  second  mode 


91 


Figure  38.  Top  section  of  block  diagram  for  second  program  mode 


Figure  39.  Bottom  section  of  block  diagram  for  second  mode 
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